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The  feasibility  instrument  was  constructed  by  modifying  existing  18-size  single-degree-of-freedom 
gyro  hardware.  Feasibility  testing  has  demonstrated  the  basic  potential  of  the  multisensor  concept  for 
avionics,  tactical  guidance,  and  navigation  applications.  Although  not  optimized  as  a multisensor,  the 
feasibility  device  performed,  as  follows,  in  a two  position,  one  hundred  hour  stability  test. 

(1)  Gyro  Bias  stability  of  0.1  meru  (0.0015  deg/hr)  and  acceleration  sensitive  drift 
of  0.3  meru/g  (0.0045  deg/hr/g). 

(2)  Accelerometer  Bias  stability  of  50  fig  and  scale  factor  stability  of  75  ppm.  The 
dynamic  and  thermal  characteristics  of  the  feasibility  instrument  are  well 
understood. 

To  obtain  similar  dynamic  performance  about  the  three  axes,  the  projected  10-size  multisensor 
utilizes  very  low  viscosity  damping  fluid.  The  active  control  is  wound  and  connected  so  that  suitable 
dynamic  performance  is  achieved,  although  the  multisensor  no  longer  functions  as  a traditional  single- 
degree-of-freedom  device.  Other  advantages  of  the  low-viscosity  fluid  are  quick  warm-up  reaction 
times,  pull-in  times  in  the  order  of  one  second,  and  rapid  instrument  fill.  With  increased  pendulosity 
and  advantages  gained  from  geometric  scaling  of  thermal  errors,  the  accelerometer  performance  should 
improve  significantly. 

In  addition  to  the  multisensor,  the  demonstrated  ability  to  sum  torque  about  a transverse  axis 
forms  the  foundation  for  a family  of  multifunction  devices  which  could  be  constructed  from  similar 
hardware.  Thus,  economic  development  and  fabrication  of  multisensors,  two-degree-of-freedom 
gyros,  two-degree-of-freedom  accelerometers,  three-degree-of-freedom  angular  rate  sensors,  and 
multisensors  is  possible. 
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SECTION  1 


INTRODUCTION 


1.1  GENERAL 

Inertial  systems  must  sense  six  quantities  (three  angular  posi- 
tion and  three  acceleration) . Current  mechanization  of  such  systems 
requires  six  single-degree-of-freedom  instruments  or  a combination  of 
two-degree-of-f reedom  and  single-degree-of-freedom  instruments  (total- 
ling four  or  five  instruments) . Multisensors  are  currently  being 
developed  at  CSDL  and  their  successful  development  will  reduce  the 
number  of  instruments  needed  for  navigation  to  three,  or  possibly  two, 
thus,  both  platform  size  and  total  system  cost  will  be  reduced.  The 
multiple  outputs  could  also  be  used  to  increase  system  reliability. 

The  multisensor  program  leads  directly  into  a family  of  multi- 
function instruments  which  will  decrease  the  costs  of  navigation 
components  and  systems.  The  multisensor  also  contributes  immediately 
to  high  accuracy  gyro  technology.  These  benefits  will  be  reviewed  after 
sinqle-deqree-of-freedom  instruments  and  tneir  extension  to  multisensors 
are  discussed. 

1.2  SINGLE-DEGREE-OF-FREEDOM  INSTRUMENTS 

A floated  single-degree-of-freedom  gyro  consists  of  a case,  a 
float,  a wheel,  and  electromagnetic  sensors  and  torquers.  The  wheel  is 
spun  at  high  angular  velocity  and  is  mounted  rigidly  to  the  float  (see 
Figure  1-1).  The  float-case  gap  is  filled  with  viscous  fluid  having  a 
density  such  that  almost  no  residual  gravity-buoyancy  force  is  exerted 
on  the  float.  Small  mismatches  in  fluid  density  and  float  unbalance  are 
countered  by  a magnetic  suspension  which  centers  the  float  with  respect 
to  the  case.  Angular  velocity  of  the  case  (hence,  torques  on  the  wheel) 
about  the  input  axis  (see  Figure  1-1)  interacts  with  the  wheel  angular 
momentum  to  rotate  the  float  about  the  output  axis.  This  angular  rota- 
tion is  sensed  by  an  electromagnetic  signal  generator,  and  is  used  in 
feedback  loops  with  a stabilized  platform  or  with  an  electromagnetic 
torquer  mounted  inside  the  instrument  case  to  maintain  the  float  at  null 
position.  The  instrument  output  is  either  the  stable-platform  position 
or  the  torque  required  to  null  the  float. 
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As  illustrated  in  Figure  1-1,  the  CSDL  multisensor  is  a single- 
degree-of-freedom  floated  gyroscope  with  pendulosity  intentionally 
added  along  the  gyro  output  axis.  Although  a cylindrical  float  is 
depicted  in  Figure  1-1,  analyses  indicate  that  noncylindrical  shapes 
merit  consideration.  The  gyro  functions  conventionally,  while  — because 
of  the  added  pendulosity  — specific  force  inputs  along  the  spin  or 
input  axes  produce  torques  about  the  input  and  spin  axes  respectively. 
These  applied  torques  are  rebalanced  by  the  float  suspension  to  main- 
tain a null  position.  The  rebalance  torques  are  measurement  outputs 
which  are  proportional  to  the  input  specific  forces. 


NOTE  THAT  EITHER  ONE  OR  BOTH  OF  THE 
ACCELEROMETER  CHANNELS  CAN  BE  USED. 


Figure  1-1.  Multisensor  concept. 


The  CSDL  multisensor  is  simple;  conceptually,  the  instrument  per- 
forms no  functions  beyond  those  required  of  a conventional  single- 
degree-of-freedom  gyro.  The  gyro  functions  as  a traditional  single- 
degree-of-f reedom  gyro,  while  mass  unbalance  has  always  been  countered 
by  suspension  forces.  The  main  difference  between  the  multisensor  and 
conventional  instruments  is  the  greater  output-axis  unbalance  and  the 
additional  measurements  which  must  be  taken  to  sense  specific  force. 
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Autonotics,  Northrop,  and  Litton  have  investigated  unrelated 
multisensor  concepts.  Autonetics  employed  an  electrostatically  suspended 
two-deqree-of-f reedom  gyro.  A tuned  rotor  suspension  with  an  essentially 
independent  gyro  and  accelerometer  was  studied  by  Litton.  The  Northrop 
device  used  quartz  crystals  for  sensing  a pressure  differential  in  liquid 
mercury  to  sense  gyro  function,  and  used  compressed  crystals  to  support 
the  instrument  and  sense  acceleration.  Teledyne,  also,  is  actively  pur- 
suing multisensor  development. 

1.3  MULTIFUNCTION  INSTRUMENTS 

Traditionally,  sinqle-degree-of-freedom  instruments  have  utilized 
the  cylinder's  axis  of  symmetry  as  the  measurement  axis.  Rather  than 
merely  demonstrating  an  added  accelerometer  function,  the  multisensor 
uses  the  transverse  axis  as  a measurement  axis.  Considerable  fabrica- 
tion and  maintenance  cost  savings  could  be  realized  from  the  following 
devices,  assembled  from  one  design  with  only  modification  to  the  float: 

(1)  Three-deqree-of-f reedom  rate  sensor  where  the  float  is 
constructed  with  no  rotating  parts  and  the  center  of  buoyancy 
nominally  coincides  with  the  center  of  gravity. 

(2)  Two-degree-of-f reedom  gyroscope,  in  which  the  float  contains 
a spinning  angular  momentum  generator. 

(3)  Two-degree-of-f reedom  accelerometer  whose  float  contains  no 
rotating  members  but  incorporates  a mass  unbalance  to  effect 
an  acceleration  sensitive  response. 

(4)  Multisensor  which  features  both  rotating  wheel  and  mass 
unbalance. 

The  multisensor  program  is  contributing  directly  to  the  develop- 
ment of  high  performance  gyroscopes.  To  realize  loop  response  required 
for  navigation,  the  multisensor  incorporates  active  suspensions  which  are 
planned  for  use  in  advanced  gyroscopes.  Multisensor  analyses  and  experi- 
ments study  radial  forces  and  motions,  phenomena  which  must  be  understood 
in  order  to  design  advanced  third  and  fourth  generation  single-degree-of- 
freedom  instruments.  In  designing  multisensors,  damping  constants  and 
their  effects  are  carefully  considered,  thereby  gaining  insight  for  other 
programs  which  seek  to  change  damping  by  decreasing  fluid  viscosity  or  by 
increasing  the  fluid  gaps. 
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This  program  combined  analyses  and  testing  of  a feasibility  instru- 
ment in  order  to  predict  the  performance  of  an  optimized  multisensor. 

The  following  were  used  as  guidelines: 

(1)  Traditional  gyro  drift  stability  less  than  one  meru 
(0.015  deg/hr) ; 

(2)  Accelerometer  stability  of  100  pg  bias  and  50  ppm  scale 
factor; 

(3)  Strapdown  application  with  maximum  case  rate  of  300  deg/s. 

The  feasibility  device  was  constructed  under  internal  funding, 
from  existing  18  IRIG  Mod  D hardware,  to  minimize  cost  and  permit  test- 
ing in  a short  time  scale.  The  feasibility  device  is  not  optimized  in 
any  sense  and  is  merely  intended  to  enable  projection  of  optimized 
configurations.  Previous  testing  was  minimal  because  of  suspected 
contamination  in  the  fluid.  Except  for  new  flex  leads,  the  removal  of 
contamination,  and  a fluid  change,  the  instrument  was  unchanged  from 
previous  builds. 

Previous  testing  and  analytic  work  resulted  in  the  following 
changes  external  to  the  instrument  prior  to  test: 

(1)  The  axial  suspension  was  made  active  to  minimize  reaction 
torques  (Paragraph  3.7); 

(2)  The  accelerometer  and  gyro  loops  were  built  to  implement  the 
desired  control  laws  and  enhance  experimental  flexibility; 

(3)  The  instrument  mounting  arrangement  was  redesigned  to  assure 
end-to-end  thermal  symmetry  (Paragraph  3.4). 

1.4  ORGANIZATION  OF  REPORT 

Section  2 describes  the  feasibility  instrument  and  test  results. 

In  Section  3,  the  tests  results  and  analyses,  which  are  included  as 
appendices,  are  combined  to  clarify  the  principal  factors  contributing 
to  multisensor  performance.  Section  4 presents  the  design  recommended 
for  a 10  size  (approximately  one  inch  case  diameter)  instrument.  The 
report  is  summarized  in  Section  5 with  recommendations  for  further 
testing  which  should  be  accomplished  before  finalizing  the  multisensor 
design. 
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SECTION  2 


FEASIBILITY  MULTISENSOR  TEST 


2.1  SUMMARY 

The  feasibility  multisensor  test  program  was  designed  as  an 
experimental  investigation  rather  than  a qualification  test,  although 
performance  guidelines  were  stated  (see  Section  1).  Understanding  the 
instrument's  behavior  must  be  emphasized.  Analyses  have  shown  that 
the  present  feasibility  instrument  is  probably  not  an  optimal  configura- 
tion. That  knowledge  coupled  with  analysis  (see  Section  4)  will  indi- 
cate improvements  and  the  potential  of  future  multisensor  designs. 

The  performance  objectives  are  1 meru  (0.015  deg/hr)  gyro  drift 
stability,  and  100  ,.g  bias  and  50  ppm  scale  factor  stability  for  the 
accelerometer  function.  The  100  hour  drift  stability  test  acceler- 
ometer uncertainties  about  the  mean  are  80  ppm  scale  factor  and  50  ng 
bias.  The  gyro  bias  drift  stability  is  0.09  meru  (0.0015  deg/hr)  and 
acceleration  sensitive  drift  stability  is  0.29  meru/g  (0.0045  deg/hr/g) . 
Stability  data  is  discussed  in  Paragraph  2.6. 

All  priority  tests  listed  in  Reference  8 were  completed  except 
for  closed  loop  frequency  response  measurements  and  noise  errors 
within  the  instrument's  bandwidth. 

Tests  which  demonstrated  the  multisensor  parameter  sensitivity 
to  changes  in  float  position,  temperature,  wheel  power,  excitation 
voltage  and  current  are  described  in  Paragraph  2.5. 

In  addition  to  parameter  stability,  a 72  hour  one  position  stability 
test  is  discussed  in  Paragraph  2.7.  Calculations  of  drift  uncertainty 
results  were  made  using  data  which  includes  trends.  After  the  trends 
are  removed,  improvements  of  0-40  percent  in  parameter  stability  are 
realized . 

A warm-up  characterization  test  was  nerformed  from  only  one  level 
of  temperature  dormancy.  The  test  results  including  instrument  tempera- 
ture are  presented  in  Paragraph  2.8. 


5 


Multisensor  cross-coupling  tests  have  been  conducted  to  evaluate 
the  effects  of  misalignment  since  the  accelerometer's  input  axis 
probably  will  not  be  on  the  reference  axes  or  coincident  with  the  gyro 
input  axis.  Also,  the  accelerometer  outputs  contain  rate  information 
(Paragraph  3.2  and  3.3).  The  data  reduction  algorithms  are  available, 
however,  since  some  difficulty  has  been  experienced  in  transferring  the 
raw  data  to  a new  computer  system,  the  results  of  this  test  are  not 
available  at  this  time. 

Tests  were  conducted  to  determine  the  cause  of  the  accelerometer's 
long  settling  time  which  was  previously  reported  (October  1976)  at 
Quarterly  Review.  The  tests  revealed  that  the  readout  active  filters 
require  up  to  30  minutes  for  the  output  to  reach  its  final  value. 

The  settle  time  is  related  to  the  initial  value  of  the  filter's  input. 
During  instrument  tumbling,  large  rate  torques  cause  the  filter  to 
saturate  or  nearly  saturate.  The  filter's  settling  characteristic  is 
attributed  to  capacitance  dielectric  which  is  somewhat  analogous  to 
magnetic  disaccommodation.  A mylar  capacitor  which  replaced  the  original 
tantulum  capacitor  eliminated  this  problem.  The  instrument  response  with 
mylar  filter  capacitor  is  given  in  Paragraph  2.8. 

2.2  INSTRUMENT  DESCRIPTION 

The  feasibility  multisensor  is  a single-deqree-of-freedom  floated 
gyroscope  with  pendulosity  intentionally  added  along  the  gyro  output 
axis . 

The  instrument  closely  resembles  the  schematic  shown  in  Figure 
1-1,  and  is  depicted  in  cutaway  form  in  Figure  2-1.  Pertinent  charac- 
teristics are  tabulated  in  Table  2-1.  The  feasibility  instrument  was 
designed  to  minimize  cost  and  to  produce  data  in  a short  time  frame; 
thus,  the  instrument  is  an  Inertial  Reference  Integrating  Gyroscope 
(18  IRIG  Mod  D)  with  modifications.  As  discussed  in  Chapter  4,  optimized 
configurations  will  probably  differ  from  the  feasibility  instrument. 
Pendulosity  was  added  by  reworking  the  balance  rings.  The  tapered 
suspension  was  replaced  by  separate  radial  and  axial  suspensions.  The 
gyro  signal  and  torque  generators  were  redesigned  to  yield  accelerom- 
eter loop  displacement  signals  (in  addition  to  performing  their  normal 
functions)  so  that  active  control  of  the  accelerometer  rebalance  could 
be  implemented.  End  housings  were  altered  so  that  separate  leads  (see 
Figure  2-2)  for  each  winding  are  available  to  increase  experimental 
flexibility.  The  float-case  gap  was  tapered  to  reduce  the  rotational 
damping  coefficients  about  the  spin  and  output  axes. 
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Figure  2-1.  Feasibility  multisensor. 

Figure  2-2  shows  the  multisensor  mounted  in  the  test  fixture 
which  was  designed  so  that  the  instrument  will  be  thermally  symmetric 
from  end  to  end  in  order  to  reduce  accelerometer  scale  factor  and  bias 
instabilities . 

2.3  BUILD  INTEGRITY 

After  discovery  or  a broken  pivot  in  the  signal  generator  end, 
the  instrument  was  refabricated  and  returned  to  test  on  24  September 
1976.  The  multisensor  was  wired  and  connected  with  passive  suspensions 
to  determine  the  acceptability  of  the  rebuilt  instrument.  All  measured 
parameters  are  within  the  specifications  of  a normal  (18  series)  gyro 
except  flotation  temperature  is  10°F  higher  than  specified  (132  F)  in 
Table  2-1.  This  flotation  discrepancy  resulted  because  the  flotation 
rings  were  removed  from  the  float  while  the  instrument  was  filled  with 
low  viscosity  fluid  during  a previous  IR&D  program.  The  measured  out- 
put axis  damping  is  reduced  by  a factor  of  2 as  a result  of  the  increased 
operating  temperature. 
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SYMFTRICAL 

MOUNTING 

FLANGE 


SG  END  HOUSING 
(NOTE  SEPARATE  LEADS 
FOR  EACH  WINDING) 


FLANGE  TEMPERATURE 
CONTROL  ANL'  MONITOR 
SENSORS 


FLANGE  HEATERS 
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L iqure  2-2.  Feasibility  multisensor  mounted  in  test  fixture. 
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TABLE  2-1. 


BUliiO  INTEGRITY  TEST  PARAMETERS. 


Axial  centering  ratio 

= 0.93 

Radial  centering  ratio 

SS  = 0.56 

SI  = 0.66 

TS  = 0.90 

TI  = 0.83 

Gyro  SG-TG  null  coincidence 

< 1.0  mr 

OA  damping  (C^,) 

= 1.14  x 10 5 dc  s 

Gyro  SG  sensitivity  with 
0 = 4 at  6 mA/leg  (I  ori) 

= 95.8  mv/mr 

OA  float  stop  angles 

+15.3  mr , -12.1  mr 

Wheel  power  (total) 

= 9.11  watts 

ASG  electrical  centers  (nulls) 
were  determined 

Float  freedom  ramp  test 
indicated  no  problems 

Flotation  along  I axis 

142. 6°F 

Flotation  along  0 axis 

144. 8°F 

When  the  accelerometer's  active  suspension  control  loons  were 
closed,  the  output  signals  were  very  noisy.  The  spin  axis  forcer 
(which  sums  torques  about  IA)  signals  exhibited  noise  level  greater 
than  the  loon  capability  of  2.5  gravities.  In  particular,  the 
forcers  on  the  float  TG  end  showed  a greater  rms  noise  level  than 
the  SG  end.  This  observation  verified  that  the  gyro  torque  rebalance 
current  interacts  magnetically  with  the  accelerometer's  signal 
generator  since  both  share  the  same  electromagnetic  circuit.  Also, 
small  signal  generator  noise  propagates  into  large  forcer  noise 
because  of  the  large  gain  required  to  achieve  wide  bandwidth  with 
high  damping  about  the  transverse  axes  (see  Section  3.6).  To 
allow  further  testing  of  the  multisensor,  the  accelerometer  band- 
width (refer  to  Appendixes  B,  C,  and  D)  was  reduced  from  22  to 
0.8  Hz  and  a small  transformer  was  added  which  couples  the  gyro 
torauer  secondary  and  the  accelerometer's  signal  generator  to 
•cancel  the  magnetic  coupling  which  occurs  within  the  instrument. 
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2.4 


DATA  PROCESSING 


The  multisensor  data  was  processed  on  a digital  computer.  After 
the  data  is  recorded  on  maqnetic  tape  using  the  Data  Acquisition  System 
(DAS),  it  is  mounted  on  the  tape  deck  of  the  PDP8  computers.  The  tape 
is  read  by  a program  which  converts  files  produced  by  the  DAS  into  files 
which  can  be  used  by  the  PDP8 . 

The  analog  signal  files  are  the  input  to  another  program  where 
they  are  transformed.  Each  raw  data  point  X is  modified  by  the  equa- 
tion Y = A X + B.  The  normal  values  of  the  analog  signals  which  have 
been  suppressed  by  passing  each  accelerometer  output  through  a Fluke 
differential  voltmeter  are  reintroduced.  Also,  each  raw  data  point  is 
scaled.  The  scaling  constant  A (Sections  3.3  and  3.9)  was  determined 
by  a calibration  test  for  each  accelerometer  axis  to  minimize  off- 
flotation-effects.  The  test  consisted  of  producing  a small  gravity 
change  to  the  accelerometer  input  axis  and  measuring  the  corresponding 
change  in  output  force  voltage  on  each  end  of  the  float.  The  constant 
is  then  calculated  by  the  ratio  of  the  two  voltages  for  each  axis. 

These  constants  reduce  the  errors  between  the  two  voltage  outputs  of 
an  accelerometer  axis  because  of  mismatched  electromagnetic  windings 
and  readout  electronics.  The  value  of  B for  a given  axis  is  the 
suppressed  voltage  for  that  axis  times  A.  The  resulting  transformed 
files  were  processed  by  the  following  options: 

(1)  Statistics  - mean,  standard  deviation,  and  coefficient  of 

vibration. 

(2)  Linear  regression  - drift  analysis. 

(3)  Plotting  - graphical  representation  of  data. 

(4)  Average  points  - filters  the  data. 

(5)  Combining  files  - Sum  of  force  signals  for  each  accelerom- 

eter. 

Four  mean  values  of  force  signals  EgI,  ETI,  Egg,  and  ETg  are 
obtained.  To  compute  scale  factor  and  bias  for  each  accelerometer, 
mean  values  are  processed  from  a two  point  test  (2  from  IA  up  and  2 
from  IA  down)  . The  values  are  an  input  to  a program  that  computes  SF 
and  bias  according  to  the  following  equations. 
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1 Accelerometers 
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where: 

EgI  = transformed  voltage  SG  end  1 axis 
Eti  = transformed  voltage  TG  end  I axis 


For  the  S accelerometers:  Substitute  9 = 0 - tt/2 ; and 

E c = transformed  voltage  SG  end  S axis  for  E„T 

OO  SI 

ETS  = transformed  voltage  TG  end  S axis  for  ETI 

Gyro  bias  and  acceleration  sensitive  drifts  are  computed  by: 


Bias  Drift  (BD) 


( E , + E 

2SF  ' down  up 


Accelerometer  Drift  = =—  (E , - E ) 

2SF  down  up 

Where  gyro  SF  is  determined  from  an  OA  vertical,  IA  horizontal 
north  to  IA  east  test 


„„  _ Enorth  Eeast 
738.9  meru 

2.5  SENSITIVITY  TESTS 

Sensitivity  testing  was  conducted  to  provide  a basis  for  compari- 
son with  stability  data  and  a better  understanding  of  the  feasibility 
instrument’ 8 behavior.  These  tests  demonstrated  the  multisensor 
sensitivity  to  the  following  conditions: 

(1)  Temperature  - uniform  and  gradient  (end  to  end). 

(2)  Wheel  power. 

(3)  SG  primary  excitation. 
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(4)  Radial  and  axial  position,  and 

(5)  Forces  primary  current. 

Sensitivity  test  results  are  listed  in  Table  2-2.  All  test  data 
is  plotted  in  Figures  2-6  through  2-35. 

The  uniform  temperature  sensitivity  was  conducted  by  varying  the 
case  temperature  controller  * j.8°F.  The  end  housing  temperature 
monitors  indicated  the  same  temperature . change  occurred  end  to  end, 
therefore,  no  gradient  was  introduced.  The  scale  factor  and  bias 
measurements  were  made  after  the  temperature  had  settled  at  least  1/2 
hour.  An  interesting  result  of  this  test  is  that  the  S-axis  SF  sensi- 
tivity is  2.3  times  greater  than  the  I-axis.  The  S-axis  bias  is  1.8 
times  larger  than  I-axis  and  of  opposite  slope.  The  test  was  performed 
twice  with  similar  results.  The  differences  between  axes  are 
attributed  to  non-ideal  temperature  distribution  and  the  voltage  scal- 
ing factors.  Further  testing  is  required. 

A thermal  gradient  test  was  performed.  A disc  type  heater  was 
cemented  to  the  SG  end  housing  to  produce  an  end  unbalance  heater  power. 

An  increase  in  SG  end  temperature  (0.072°F)  caused  the  case  tempera- 
ture controller  to  respond  with  a decrease  and  a reduced  temperature 
in  the  TG  end  -0.02°F.  The  new  result  is  a nonlinear  temperature 
gradient  along  the  instrument's  OA  as  shown  in  Figure  2-3.  Note  that 
the  end  mounted  thermistors  indicated  a 1.5UF  cooler  than  normal  operating 
temperature  than  the  center  of  the  instrument . Because  the  OA  thermal 
gradient  is  not  linear,  the  sensitivities  for  this  test,  listed  in 
Table  2-2,  should  be  interpreted  as  approximate  values  only. 

Wheel  power  sensitivity  was  determined  by  varying  the  excitation 
+ 5 percent  from  a normal  of  9.1  watts.  Sufficient  time  was  allowed 
for  temperature  settling  before  scale  factor  and  bias  measurements 
were  made.  A 5 percent  change  in  wheel  power  causes  a 0.02°F  change 
in  instrument  case  temperature.  The  temperature  controller  responds 
to  decreasing  wheel  power  by  increasing  heater  power.  The  S-axis 
scale  factor  sensitivity  is  4.5  times  larger  than  along  the  I-axis. 

The  difference  between  axes  is  attributed  to  uncertainty  in  tempera- 
ture distribution  and/or  in  the  voltage  scaling  factors. 

The  accelerometer  scale  factor  sensitivities  to  SG  primary  exci- 
tation are  small  and  nearly  identical  in  both  axes.  The  S-axis  bias 
sensitivity  is  20  times  larger  than  the  I-axis  sensitivity.  Similarly, 
the  reaction  torque  attributed  to  the  forcer  primary  is  20  times 
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TABLE  2-2.  FEASIBILITY  INSTRUMENT  SENSITIVITY  TEST  RESULTS 


larger  in  the  S-axis  than  in  the  I-axis.  These  three  observations  are 
explained  by  misalignments  between  the  zero  torce  positions  of  the 
signal  generator,  the  forcer  primary  zero  force  position,  and  the 
signal  generator  zero  displacement  signal  position. 


+0A 


DISTRIBUTION 


Figure  2-3.  Multisensor  normal  and  temperature  gradient  distribution. 


In  preliminary  air  bearing  tests  at  the  component  level,  the 
accelerometer  signal  generators  were  tuned  with  positive  elastic 
restraint  to  balance  the  negative  restraint  of  the  constant  magnitude 
square  wave  forcer  primaries.  Elastic  restraint  is  defined  as  a posi- 
tion sensitive  force.  This  technique's  success  is  reflected  in  the  I- 
axis  bias  sensitivity  to  I-axis  radial  position  (250  pg/pin,  not  very  qood 
since  this  is  near  the  restraint  of  forcer  primary  alone)  and  similarly  for 
the  S-axis  (10  pg/min,  an  acceptable  value)  (Appendix  D elaborates). 

To  account  for  assymmetries  between  components  and  other  sources  of  elastic 
restraint  such  as  the  axial  suspension,  the  elastic  restraint  tuning 
should  be  done  in  the  actual  instrument  as  will  be  discussed  shortly. 

Return  to  the  alignment  of  the  various  null  conditions.  Because 
of  the  integrator  in  the  accelerometer  loops,  the  float  position  with 


14 


respect  to  the  case  is  determined  when  the  output  of  the  signal  genera- 
tor is  null.  With  current  transformers,  this  corresponds  to  equal 
currents  on  opposite  magnetic  cores.  The  signal  generator  is  tuned  as 
a passive  suspension  so  that  the  zero  force  point  which  depends  on  the 
capacitors  may  differ  from  the  equal  current  position.  The  sensitivity 
to  signal  generator  excitation  is  proportional  to  the  displacement 
between  the  two  centers.  If  the  centers  do  not  coincide,  changing  the 
excitation  causes  the  signal  generator  to  exert  a force.  The  controller 
then  alters  the  forcer  secondary  current  so  that  the  equal  current  con- 
dition is  met  and  the  sum  of  forces  exerted  on  the  float  is  zero. 
Similarly,  a displacement  between  the  SG  zero  current  position  and 
the  forcer  primary  null  force  point  results  in  the  reaction  torque 
attributed  to  the  forcer  primary. 

A further  explanation  of  the  situation  and  estimate  of  the  degree 
of  misalignment  are  offered  in  Figure  2-4  which  models  the  forcer  and 
SG  at  one  end.  The  letters  refer  to  the  following: 

(1)  Line  A is  the  force  displacement  characteristics  of  the 
forcer  primary. 

(2)  Line  B is  the  force-displacement  characteristic  of  the 
signal  generator. 

(3)  Point  C is  the  zero  force  point  for  the  signal  generator. 

(4)  Point  D is  the  zero  force  for  the  forcer  primary. 

(5)  Axis  E is  the  equal  current  null  position  of  the  signal 
generator.  Because  of  the  integral  control  loop,  steady- 
state  operation  will  result  on  this  axis. 

(6)  Force  F is  the  reaction  of  the  forcer  primary. 

(7)  Force  G is  the  reaction  of  the  signal  generator. 

(8)  Force  H is  the  force  exerted  by  the  forcer  secondary  so 
that  the  sum  of  forcer  (H  plus  F)  and  signal  generator 
forces  (G)  is  zero  (with  no  external  forces  operating) . 

(9)  Slope  I is  the  elastic  restraint  of  the  forcer. 

(10)  Slope  J is  the  elastic  restraint  of  the  signal  generator. 

The  situation  in  the  actual  instrument  is  complicated  by  the 
presence  of  two  ends;  however,  the  basic  argument  remains  valid.  For 
the  forcer  restraint  of  1.6  * 105ugand  measured  elastic  restraint 
(air  bearing  tests)  of  2 « 10 5 dyn-cm/cm,  the  displacement  between 
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signal  generator  null  and  forcer  null  is  260  gin.  a fiqure  certainly 
nossible  with  manufacturing  tolerances. 


FORCE 


POSITION 


Figure  2-4.  Performance  of  SG  and  forcer. 

(see  text  for  nomenclature) 

In  future  testing,  bias  and  elastic  restraint  tuning  will  be 
performed  in  the  instrument.  For  the  present  device  where  the  signal 
generator  is  connected  in  series  parallel  connection,  elastic  restraint 
can  be  adjusted  by  increasing  the  capacitance  on  both  legs  of  a given 
axis  simultaneously.  Bias  can  be  adjusted  by  either  of  the  following 
techniques : 

(1)  Change  the  SG  null  by  adding  a dummy  signal; 

(2)  Push-pull  the  SG  capacitors;  and/or 

(3)  Add  shunt  resistors  parallel  to  the  legs  of  the  forcer 
primary . 

Since  future  designs  will  arrange  the  SG  legs  in  bridge  circuits, 
option  2 can  be  implemented  by  adding  shunting  capacitors  parallel  to 
the  inductances. 

The  instrument's  float  position  sensitivity  was  determined  by 
dummy  directing  the  active  suspension  loops  along  each  of  the  S-,  I-, 
O-axes.  Sufficient  time  was  allowed  for  transients  to  settle  before 
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taking  data.  Axial  and  cross-axis  displacement  scale  factor  sensitivi- 
ties are  less  than  10  ppm/tiin.  I-axis  and  S-axis  scale  factor  as  a 
function  of  position  along  their  respective  axes  display  a distinct 
nonlinear  characteristic  particularly  with  I-axis  position  as  shown  in 
Figure  2-28.  The  test  was  repeated  with  similar  results;  however, 
more  data  points  will  be  required  to  fully  characterize  the  instrument's 
behavior.  The  nonlinearity  in  scale  factor  could  be  caused  by  non- 
linearity in  the  force  position  relationships  of  magnetic  components. 


2.6  STABILITY  TESTS 


The  multisensor  was  tested  for  a period  of  100  hours  concentrating 
on  the  accelerometer  scale  factor  and  bias  stability.  The  instrument 
was  positioned  OA  horizontal  east,  parallel  to  the  table,  with  IA 
up  and  down  perpendicular  to  the  earth’s  axis  to  minimize  w0A  coup- 
ling (see  Figure  2-5  and  Paragraph  3.2).  Ac  each  position,  sufficient 
time  was  allowed  for  the  forcer  signals  to  settle  before  accumulating 
data.  Ten-minute  segments  of  one  second  sampled  data  is  recorded  on  a 
data  acquisition  system  magnetic  tap  for  computer  analysis.  Additional 
measurements  recorded  during  the  stability  tests  monitored  the  following: 


(1)  Instrument  end  housing  and  case  temperature; 
(21  Changes  in  wheel  power; 

(3)  Axial  force,  and 

(4)  Primary  excitation. 


Figure  2-b.  Multisensor  position  for  scale  factor 
and  bias  calibration  tests. 
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The  stability  test  results  and  performance  objectives  are 
listed  in  Table  2-3.  All  data  is  plotted  as  shown  in  Figures  2-36 
through  2-44.  These  results  indicate  that  excellent  gyro  performance 
has  been  obtained.  The  bias  drift  standard  deviation  exceeds  the  goal 
(0.015  deg/hr)  by  at  least  a factor  of  ten  (0.00135  deg/hr).  This 
performance  is  attributed  to  the  independent  active  control  of  the 
radial  and  axial  suspensions  despite  the  large  OA  pendulosity  (by  gyro 
standards) . 

The  accelerometer  performance  data  is  within  the  goals  except  for 
the  marginal  SF  (62.3  ppm)  uncertainty  of  the  S-axis.  Wheel  power  and 
SG  primary  excitation  stability  measurements  do  not  contribute  signifi- 
cantly to  the  accelerometer  uncertainties.  Wheel  power  stability  is 
<0.002  watts  contributing  a maximum  uncertainty  of  8 ppm  scale  factor 
based  on  independent  sensitivity  measurements  listed  in  Table  2-2. 

The  SG  primary  excitation  stability  was  measured  at  0.045  percent  but 
due  to  equipment  difficulties  at  the  time  the  source  was  not  monitored, 
during  the  stability  test.  The  contribution  of  this  source  to  scale 
factor  uncertainty  is  1.5  ppm.  The  S-axis  has  a major  bias  sensitivity 
(65  vig/0.045  percent)  to  primary  excitation  while  the  I-axis  is  20 
times  smaller  (3.2  ug/0.045  percent).  The  bias  standard  deviations  are 
nearly  the  same  for  both  axes  during  stability  testing;  therefore,  it 
is  reasonable  to  assume  then  that  this  excitation  did  not  change 
significantly  during  the  test. 

Multisensor  temperature  variation  is  a major  source  of  scale 
factor  instability  according  to  analysis  and  sensitivity  measurements, 
(see  Paragraphs  2.5  and  3.4).  Based  on  correlation  of  thermal  measure- 
ments with  accelerometer  data,  temperature  was  the  major  contributor 
to  the  stability  test  uncertainties.  A careful  examination  of  test 
data  plotted  in  Figures  2-36,  2-37  and  2-42  shows  a definite  agreement 
in  scale  factor  and  temperature  characteristics.  In  addition,  the 
-0.009  F temperature  drift  will  produce  approximately  10  ppm  drift  in 
scale  factor.  Also,  scale  factor  variations  about  the  trend  line 
have  relative  magnitudes  that  agree  well  with  temperature  oscillations. 

The  temperature  plot  of  Figure  2-42  was  produced  from  strio 
chart  recordings  with  the  instrument  IA-SA  positioned  down.  This  plot 
is  also  representative  of  the  other  three  thermistor  characteristics 
which  are  not  presented.  Instrument  case  temperature  is  shown  plotted 
in  Figure  2-43.  Despite  some  difficulty  with  a noisy  Sanborn  strip 
chart  recording , particularly  during  the  first  half  of  the  test,  the 
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70.5  raeru/g  (unbalance  along  SA) 


case  temperature  oscillations  jyiee  with  the  end  housing  measurements 
from  72  hours  to  100  nours  ol  the  test.  The  case  temperature  positive 
trend  is  due  to  the  temperature  controller  reaction  to  a decreasing 
trend  in  ambient  temperature  as  indicated  by  base  and  end  housing 
measurements . 

2.7  ONI'- POSIT  I ON  STABILITY  TEST 

The  multisensor  was  operated  with  output  axis  horizontal  and 
east  with  the  input  axis  up  and  perpendicular  to  earth  axis  (42  degrees 
from  vertical)  for  72  consecutive  hours.  A commercial  tilt  meter 
mounted  on  the  test  stand  measured  tilt  angles  in  the  north-south  and 
east-west  planes.  The  IA-SA  plane  lies  in  the  north-south  plane. 

Data  was  recorded  on  the  data  acquisition  system.  The  four 
radial  forcers,  the  gyro  torquer,  four  case  thermistors,  wheel  power 
and  the  platform  tilt  were  monitored.  Sampled  at  ten  second  intervals 
and  averaged  over  two  minutes,  the  data  versus  time  is  plotted  with 
statistical  information  in  Figures  2-45  through  2-51.  Since  the  four 
thermistors  were  nearly  identical,  only  one  temperature  plot  is 
included. 

Power  spectral  densities  are  plotted  in  Figures  2-52  through 
2-54.  Aliasing  effects  were  assumed  small  for  the  following  reasons: 

(1)  Data  averaging  acted  as  a low  pass  filter; 

(2)  The  accelerometer  channels  were  filtered  at  0.02  Hz; 

(3)  Analog  strip  chart  recordings  did  not  reveal  rising  spectral 
characteristics  beyond  the  Nyquist  frequency. 

Leakage  (the  sampling  of  a nonintegral  number  of  cycles)  has  affected 
the  spectral  plots  significantly.  The  tilt  versus  time  plots  clearly 
demonstrate  spectral  peaks  which  leakage  spreads  as  a i/f2. 

The  accelerometer  standard  deviations  about  the  trend  are  16.5  ppm 
for  the  I-axis  and  26.3  ppm  for  the  S-axis,  figures  which  are  two  and 
three  times  better  than  the  scale  factor  and  bias  measurements  of  the 
two  position  tests  described  in  Paragraph  2.6.  The  S-axis  accelerometer 
drifted  -33  ppm  during  the  72  hours  while  the  I-axis  drifted  -169  ppm. 

The  drift  of  the  two  axes  is  similar  if  the  step  at  55  hours  in  the  I- 
axis  data  is  omitted.  This  shift  could  be  caused  by  a capacitor  shift 
since  this  axis  has  a large  position  sensitivity  (250  pg/pin) . Appendix  D 
describes  this  phenomenon  more  completely.  With  the  sensitivities 
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calculated  for  averaqe  temperature  changes,  the  trend  of  0.0023  F/hr 
could  account  for  -13.8  and  -8  ppm  in  the  S-  and  1-axes  respectively. 

Both  accelerometer  channels  exhibit  low  frequency  oscillation  (less 
than  two  in  the  72  hours)  which  leak  to  form  the  accelerometer  snectrals, 

1/f2  slope  at  low  frequencies.  The  temperature  spectrum  between  10~4  Hz 
and  4 x 10  Hz  is  approximately  white  noise  of  amplitude  6 x 10-  F /Hz. 
Multiplying  the  temperature  spectrum  by  the  I-axis  average  temperature  sen- 
sitivity and  using  5.9  volts  as  one  gravity,  one  obtains  an  accelerometer 
- 7 2 

spectrum  of  2 x 10  volt  /Hz  which  is  the  value  of  the  I-accelerometer  soec- 

- 3 - 3 

trum  between  2 x 10  and  4 x 10  Hz.  Performing  a similar  calculation  on 

— 6 2 

the  S-axis  results  in  2 x 10  volt  /Hz,  a value  which  is  an  order  of  mag- 
nitude less  than  the  horizontal  portion  of  the  spectrum. 

Tne  wheel-power  spectrum  and  the  accelerometer  sensitivities 
also  contribute  to  the  accelerometer  errors.  Fitting  the  wheel 
spectrum  with  white  noise  of  0.0015  watt2/Hz,  accelerometer  spectra 
attributed  to  wheel-power  variation  are  6 x 10“ 7 volt2/Hz  for  the 
S-accelerometer  and  4 x 10-e  volt2/Hz  (see  Figure  2-58).  Both  figures 
are  an  order  of  magnitude  less  than  the  plotted  spectra.  Consideration 
of  this  sensitivity  is  not  entirely  valid  since  wheel-power  shifts 
affect  temperature  which  has  been  considered  separately.  However, 
gyro  testing  has  revealed  coupling  between  flex  leads  and  motor  wind- 
ings which  do  not  depend  on  temperature  so  that  the  wheel-power 
correlation  is  included.  Although  the  precise  cause  of  the  S-accelerom- 
eter spectrum  between  10“ 5 Hz  and  4 x 10~3  is  unclear,  the  following 
comments  are  offered: 

(1)  The  spin  axis  contains  information  about  case  rates  about 
the  output  axis. 

(2)  The  coupling  of  the  accelerometer  signal  generator  and  the 
gyro  is  probably  not  the  source.  With  white  torque  noise, 
the  erroneous  signal  in  the  accelerometer  signal  generator 
would  also  be  white  noise  because  of  the  inductive  coupling. 

From  Appendix  D,  the  torque  rebalance  noise  at  low  frequen- 
cies would  rise  proportional  to  frequency  squared  which 
disagrees  with  the  plotted  spectrum. 

(3)  Because  of  uncertainties  in  measuring  sensitivities  and 
temperature  profiles,  estimated  sensitivities  could  be  in 
error  by  a factor  of  two.  An  increase  in  sensitivity  of 
scale  factor  to  temperature  level  would  explain  the  flat 
portion  of  the  S curve. 
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The  time  data  and  spectrum  for  tilt  indicate  that  tilt  has  only 
a small  effect  on  stability  data.  Figures  2-55  and  2-56  show  an  rms 
tilt  of  2.7  * iO"6  rads  at  one  cycle  per  day  which  represents  an  accel- 
erometer error  of  only  2 ppm  when  the  axes  are  45°  to  the  vertical. 

To  summarize  the  accelerometer  data,  temperature  stability  impacts 
accelerometer  stability  greatly  although  other  factors  are  present. 

The  sensitivities  to  radial  positions  and  the  drift  data  and  low  fre- 
quency spectra  strongly  indicate  another  mechanism  such  as  capacitor 
shift.  Further  study  is  recommended. 

The  gyro  torque  spectrum,  Figure  2-54,  approximately  agrees  with 

data  compiled  from  conventional  18  IRIG  and  TGG.  The  cited  references 

attribute  the  noise  to  internal  disturbances  such  as  temperature 

gradients  across  the  float,  turbulence  in  the  air  bearing,  and  signal 

generator  noise  and  to  environmental  input.  As  indicated  in  the 

temperature  spectrum,  several  areas  of  the  temperature  spectrum  closely 

match  those  of  the  gyro;  however  attempts  to  correlate  the  two  through 

the  gyro  aca’e  factor  sensitivity  to  average  temperature  result  in  gyro 

torque  spectra  of  3 * 10“"  meru2/Hz  compared  to  the  plotted  spectra  of 
2 

0.5  meru  /Hz.  Temperature  gradient  along  a radial  axis,  a principle 
driver  of  gyro  error  (Reference  7),  was  not  monitored. 

2.8  REACTION/WARM-UP  CHARACTERISTIC 

The  purpose  of  this  test  was  to  demonstrate: 

(1)  The  thermal  time  constant  of  the  instrument,  and 

(2)  Pull-in  characteristics. 

The  reaction  time  of  the  multiBensor  was  tested  by  removing  the 
instrument  heater  power,  decreasing  the  operating  temperature  1.5°F. 

The  temperature  controller  causes  a 2.0°F  overshoot  during  warm-up  as 
shown  in  Figure  2-59.  The  test  was  performed  with  the  input  axis 
vertical  and  the  spin  axis  horizontal.  A decrease  in  measurement 
resolution  was  required  to  maintain  on-scale  recordings  of  instrument 
response.  Because  of  reduced  resolution  capability,  baseline  measure- 
ments were  conducted  for  comparison  purposes  prior  to  and  following 
the  transient  test.  The  multisensor  outputs  were  stable  within  five 
minutes  after  the  temperature  settled.  I-axis  forcer  output  plots 
(gravity  sensitive)  are  shown  in  Figures  2-60  and  2-61,  with  the 
resultant  summed  (note  concellation  by  summation)  outputs  in  Figure 
2-62.  The  S-axis  (insensitive  to  gravity)  outputs  are  plotted  in 
Figures  2-63,  2-64  and  2-65.  Multisensor  temperature  presented  in 
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Figure  2-59  is  monitored  by  one  of  the  four  thermistors  mounted  on  the 
end  housings.  The  other  thermistor  responses  are  the  same. 

The  gyro  response  from  this  test  is  not  presented  here,  because 
of  large  gyro  output  voltage  due  to  the  vertical  component  of  earth 
rate  (673.8  meru)  about  the  gyro  input  axis  and  the  DAS  automatic  rang- 
ing limited  resolution  capability  of  tl.u  millivolt  (0.5  meru).  Volt- 
age suppression  was  not  used  on  the  gyro  output  prior  to  data  acquisi- 
tion. Based  on  data  from  a trial  test  (gyro  IA  perpendicular  to  EA) , 
the  gyro  output  is  stable  within  five  minutes  after  the  temperature 
has  settled.  The  gyro  response  characteristic  durinq  warm-up  from 
this  test  is  plotted  in  Figure  2-66. 

The  multisensor's  reaction  to  axial  float  motion  was  measured  as 
the  float  returned  to  its  center  position  from  a 120  microinch  TG  dis- 
placement during  sensitivity  testing.  The  SG  end  S-axis  force  output 
and  axial  suspension  current  did  not  record  properly  on  the  DAS  during 
this  test.  Strip  chart  recording  of  the  S-axis  accelerometer  shows 
similar  results  (not  shown)  as  the  I-axis  outputs  plotted  in  Figures 
2-67  and  2-68.  Strip  chart  recordings  indicate  a 50  minute  axial 
suspension  settling  time.  The  accelerometers  and  gyro  are  stable  35 
minutes  after  the  suspension  is  turned  on.  The  gyro  response  is  plotted 
in  Figure  2-69. 

Temperature  measurements  made  on  the  SG  and  TG  end  housings  are 
shown  in  Figures  2-70  and  2-71.  Since  the  float  was  initially  in  the 
TG  end,  the  active  suspension  applied  pull-in  current  to  the  SG  end 
winding  causing  a 0.075°F  increase  in  end  temperature.  The  tempera- 
ture controller  responded  by  a decrease  in  power,  cooling  the  Tg  end  by 
0.045  F.  Approximately  16  or  17  minutes  was  required  for  the  tem- 
perature to  stabilize  after  turn-on.  The  accelerometer's  long  settling 
characteristic  is  caused  by  large  reaction  forces  due  to  the  axial  sus- 
pension's stator/rotor  relative  cocking  and/or  uncovering.  The  axial 
suspension's  large  reaction  forces  were  discovered  during  a previous 
multisensor  IR&D  program. 

Section  4 introduces  design  modification  which  will  speed  pull-in 
time  considerably. 
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Figure  2-8.  I-Axis  acceleration  bias  temperature  sensitivity 
(uniform  temperature  change)  . 


Fiqure  2-9.  S-Axis  acceleration  bias  temperature  sensitivity 
(uniform  temperature  change) . 
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Figure  2-10.  Gyro  bias  drift  temperature  sensitivity 
(uniform  temperature  change) . 
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Figure  2-11.  Gyro  acceleration  drift  temperature  sensitivity 
(uniform  temperature  change) . 
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Figure  2-12.  I-Axis  acceleration  scale  factor 
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Figure  2-13.  S-Axis  acceleration  scale  factor 

wheel  power  sensitivity  (45  percent). 


27 


I AXIS  BIAS  (Ml 


Figure  2-20.  I-Axis  accelerat ion  bias  sensitivity  to 
SG  primary  excitation  (9.6  kHz). 


Figure  2-21.  S-Axis  acceleration  bias  sensitivity  to 
SG  primary  excitation  (9.6  kHz). 
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Figure  2-24.  I-Axis  scale  factor  sensitivity  to  axial  position. 


Figure  2-25. 


S-Axis  scale  factor  sensitivity  to  axial 


position. 
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2-30.  S-Axis  acceleration  scale  factor  sensitivity  to 
S-Axis  displacement. 
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Figure  2-31.  S-Axis  scale  factor  sensitivity  to  I-Axis  displacrment 
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Figure  2-32. 
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Figure  2-33.  I-Axis  acceleration  bias  sensitivity  to 
S-Axis  displacement. 
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Figure  2-38.  I-Axis  acceleration  bias  stability  10/22  - 10/26/7? 
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Figure  2-41.  Gyro  AD  stability  10/22  - 10/26/76 
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Figure  2-42.  SG  end  housing  temperature  10/22  - 10/26/76. 
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Figure  2-44.  Base  temperature  10/22  - 10/26/76. 
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Figure  2-45.  72-Hour,  I-Axis  acceleration  drift  (ia  I EA) . 
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igure  2-4G.  72-Hour,  S-Axis  acceleration  drift  (SA  jj  EA) . 
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Figure  2-49.  72-Hour,  table  tilt,  north-south  plane 


Figure  2-50.  72-Hour,  table  tilt,  east-west  plane. 


421  POINTS  OBSERVATIONS  REGRESSION:  Y-A+BX 


-5  -4  -3  7 

FREQUENCY  <H*| 

Fiqure  2-52.  I-Axis  accelerometer  power  spectrum  (trend  removed). 
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Fiqure  2-55.  Power  spectrum  ot  Temperature  (includes  trend). 
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Figure  2-58.  Power  spectrum  of  wheel  power 
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Figure  2-62.  I-Axis  acceleration  temperature  reaction  (IA  vertical)  (12/13/76) 
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Figure  2-63.  SS  torcer  temperature  reaction  (SA  horizontal)  (12/13/76) . 
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Figure  2-67.  I-Axis  acceleration,  SG  end  forcer  (SI)  output  during  axial  position. 
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Figure  2-70.  SG  end  housing  temperature  during  axial  pull-in 
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SECTION  3 


DESIGN  CONSIDERATIONS 


* 


t 

i 

I 


3.1  GENERAL 

To  complement  the  description  of  the  test  program,  design  con- 
siderations unique  to  the  multisensor  are  introduced.  These  observa- 
tions dictate  the  design  of  Section  4 from  test  and  analysis.  The 
following  areas  have  been  identified  as  design  topics  which  are  particu- 
larly pertinent  to  the  multisensor  or  other  multifunction  devices: 

(1)  Equations  of  motion.  The  motion  equations  are  introduced 
in  a form  suitable  for  multisensor  application.  The  equa- 
tions introduce  topics  which  are  discussed  in  additional 
detail  throughout  this  chapter. 

(2)  Magnetic  interactions  between  axes  and  location  of  axes. 

(3)  Thermally  induced  error. 

(4)  Accelerometer  errors  introduced  by  wheel  anisoelasticity 
and  hunt. 

(5)  Accelerometer  resolution  considering  forcer  elastic  restraint 
and  signal  generator  noise  and  null  shift. 

(6)  Reaction  torques  caused  by  axial  suspension. 

(7)  Adjustment  of  damping  constants  - the  interaction  between 
input-  arid  output-axis  loops  coupled  through  wheel  angular 
momentum. 

(8)  Control  of  rotational  and  translational  modes. 

(9)  Case  rates  of  300  deg/s. 

3.2  EQUATIONS  OF  MOTION 

The  equations  of  motion  with  underlying  assumptions  are  derived 
in  Appendix  A and  are  repeated  here  for  convenience.  The  motion  equa- 
tions for  translational  motion  are: 


+ CTXS  ” (B"M)as  + Faus  (S)  + Fadd  (S) 


(3-1) 


mXI  + CtXi  = (B-M) aj  ♦ Fsus  (I)  + Fa<Jd  (J) 


mX_.  + C.VX^  = (B-M)a_  + F + F .. 

0 AX  O 0 sus  (O)  add  (0) 


(3-2) 

(3-3) 


For  rotation,  the  equations  of  motion  are: 


rs0s  + CR  9s  aIP  a0  UI  + MSUS  (S)  + Mg2  (S)  + Madd  (S) 


+ Mbuoy  (S,  “ H " TS  “s  + (IrV  < W (W 


(3-4) 


II9I  + CR  9I  * H0O  = “a8P  + a~  U~  + M 


'<">  US  ' “bus  (I)  + Mg2  (I)  + Madd  (I) 


+ “buoy  (I)  - H “o  “ JI  “l  + (W  {WS+0S)  (U0+O0) 


(3-5) 


I090  + C090  H0  I ~ aS  UI  " aI  US  + Msus  (O)  ~ Mg2  (O)  + Madd  (0) 

+ Mbuoy  (O)  + H “i  - J0  “0  + ‘W  (VV  (“l+0I)  <3~6) 


where  the  subscripts  S,  I,  and  0 refer  to  the  spin,  input  and  output 
axes  respectively  and  the  superscript  dot  (•)  implies  differentiation 
with  respect  to  time. 


6 - float  angular  position  measured  with  respect  to  use 

id  = angular  rate  inputs  applied  to  instrument  case 

a * acceleration  input  applied  to  instrument  case 

B = float  volume  times  fluid  density  which  is  temperature 

dependent 

damping  constant  along  output  axis 
damping  constant  about  output  axis 
damping  constant  about  radial  axis 
deunping  constant  along  radial  axis 


'AX 
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Fadd 


translational  forces  which  are  not  accounted  in  other 
terms 


F 

sus 


forces  exerted  by  magnetic  suspensions 


H = wheel  angular  momentum  where  the  spin  velocity  is  measured 
with  respect  to  the  float 

I = float  moments  of  inertia  including  wheel 

m = mass  of  float  including  wheel 


Madd  = torcIues  which  are  not  included  in  other  terms 


Mbuoy 


convection  torques  which  are  function  of  temperature 
distribution.  Since  these  torques  are  acceleration 
sensitive,  they  can  be  included  in  pendulosity  and 
unbalance  terms 


Mg2  = torques  sensitive  to  acceleration  squared 

P = pendulosity  along  output  axis 


Mgus  = torques  exerted  by  the  magnetic  suspensions  and  torquers 
U = mass  unbalance 


X = translational  position  of  the  nominal  float  center  with 

respect  to  the  case 


Many  terms  in  the  above  equations  are  discussed  in  other  portions 
of  this  chapter  and  in  the  appendices;  however  several  comments  are 
offered.  As  in  single-degree-of-freedom  devices  the  output  axis  is 
primarily  driven  by  H Wj.  the  product  of  the  wheel  momentum  and  the  case 
rate  about  the  input  axis.  Because  of  the  symmetry  of  the  motion 
equations,  the  torque  about  the  input  axis  includes  H u)Q,  the  case 
rate  about  the  output  axis.  The  torques  about  the  input  and  spin  axes 
contain  terms  which  multiply  the  pendulosity  along  the  output  axis  by 
the  accelerations  along  the  spin  and  input  axes  respectively.  The 
pendulosity  which  is  the  moment  generated  by  the  noncoincidence  of  the 
float's  centers  of  gravity  and  buoyancy  drives  the  accelerometers. 

Thus,  before  cross-couplings  are  included  (paragraph  3.3),  the  following 
observations  are  offered. 

(1)  The  output  axis  maintains  the  gyro  function; 

(2)  The  spin  axis  functions  as  an  accelerometer; 

(3)  The  input  axis  contains  both  accelerometer  and  gyro 
information. 
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Although  three  information  channels  are  available,  the  input-axis  channel 
contains  both  accelerometer  and  gyro  information  so  that  information 
from  another  sensor  is  required  to  separate  the  information.  The  system 
designer  has  the  option  of  using  three  multisensors  to  obtain  the  six 
measurements  typically  used  in  navigation  and  will  then  have  three  extra 
channels  containing  redundant  information.  With  six  variables  and  six 
channels  of  information,  theoretically,  two  multisensors  can  mechanize 
an  entire  platform. 

The  rotational  motion  equation  about  the  input  and  output  axes 
are  coupled  through  Hw  terms.  This  coupling  is  discussed  in  Section 
3.0  and  in  Appendix  C and  will  be  exploited  in  the  designs  of  Section  4. 
The  spin  axis  equation  is  essentially  uncoupled  from  the  other  rota- 
tional axes.  Implementation  of  the  magnetic  suspensions  will  couple 
the  translational  motion  along  the  spin  axes  and  the  rotation  about 
the  input  axis  and  for  translational  along  input  and  rotation  about 
spin,  a phenomena  which  is  discussed  in  Section  3.9  and  Appendix  B. 
Coupling  introduced  by  misalignments  and  magnetic  coupling  are  st”died 
in  Section  3.3.  Terms  which  include  angular  rate  products  will  be 
important  for  strapdown  applications  and  are  discussed  in  Section  3.10. 

3.3  MAGNETIC  INTERACTIONS  AND  LOCATION  OF  AXES 

The  previously  derived  equations  of  motion  assumed  that  the 
reference  axes  were  aligned  with  the  physical  spin,  input,  and  output 
axes  of  the  float.  This  section  introduces  the  consequences  of  refer- 
ence axes,  misalignment,  and  cross-couplings  of  magnetic  suspension 
and  torquers.  These  phenomena  determine  the  actual  axes'  location 
with  respect  to  the  reference  axes  and  determine  the  compensation 
required  for  navigation. 

The  reference  axes  are  lines  scribed  on  the  instrument's  case; 
float  axes  are  defined  as  the  physical  axes  of  the  float.  Without  arguing 
about  the  input  and  output  axes,  the  float  spin  axis  would  be  the  line 
about  which  the  wheel  spins.  As  derived  in  Appendix  F,  the  instrument’s 
physical  axes  are  defined  by  the  position  of  the  float  physical  axes 
and  by  magnetic  misalignments  and  cross-couplings,  mass  unbalancer, 
convection  torques,  and  buoyant  forces. 

To  successfully  navigate,  the  actual,  not  ideal,  content  of  a 
torque  channel  must  be  determined.  Even  with  perfect  components  assumed 
the  instrument  axes  would  be  misaligned  with  the  navigation  coordinates 
because  of  machining  tolerance  inherent  in  manufacture;  thus,  compensation 
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would  be  required.  For  the  ensuing  discussion  the  case  axes  are  assumed 
aligned  with  the  navigation  axes. 

In  Appendix  F,  a static  analysis  of  the  multisensor  is  conducted 
where  the  float  axes  are  assumed  misaligned  with  respect  to  the  case 
reference  axes  and  linear  magnetic  misalignments  and/or  cross-couplings 
are  modeled.  The  resulting  signals  which  would  be  ted  into  the  naviga- 
tion computer  are  summarized  in  Equations  (F-14)  and  (F-15) . These 
results  are  used  to  model  the  axes  location  tests  described  in  Section  2. 
For  earth  rate  testing,  the  product  ot  rate  terms  may  be  neglected. 

Based  on  physical  arguments  leading  to  (F-14)  and  (F-15) , the 
performance  equation  for  the  reference  gyro  output  axes  may  be  written 
as 


W 


IND 


+ wo“o  ♦ Ks“s  * 

df 

DI 

ai  + Ds 

+ Doo  *o!  * Dss 

"s'  * DI0  ai 

ao 

+ °IS  a 

(3-7) 


where: 

bp  is  a term  insensitive  to  acceleration  or  rate  (bias) 

DI ' DS'  D0  are  coe^^lcients  acceleration  sensitive  terms 

D00'  DII'  DSS'  UIO'  DIS'  Dns  are  coefficients  terms  sensitive 

to  acceleration  squared 

WI , w . Ws  are  terms  proportional  to  case  rates 
a indicates  acceleration 

indicates  angular  rate 
u>  indicates  rate  along  reference  axes 

The  subscripts  S,  I,  and  0 indicate  vectors  along  the  spin,  input, 
and  output  case  reference  axes. 

The  acceleration  terms  are  called  cross-coupling.  Ideally, 
only  the  coefficient  Wj  should  be  nonzero;  however,  the  coefficients 
Wg  and  WQ  can  be  called  misalignment  or  cross-coupling.  Assuming  the 
Wg/Wj  and  wq/wj  are  much  less  than  one,  rotate  the  reference  axes 
-W-./W  about  the  spin  axis  and  -W  /W  about  the  input  axis;  thus,  the 
input  axis  will  be  defined  in  the  new  coordinates  so  that  only  rates 
about  the  new  input  axis  are  sensed. 
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Similar  expressions  result  for  the  other  channels  where  an 
accelerometer  can  be  defined  as  aligned  when  the  linear  cross-accelera- 
tion terms  are  zero.  Since  the  reference  axes  can  only  be  rotated 
about  three  axes,  and  alignment  of  each  of  the  three  axes  requires  two 
rotations,  all  the  axes  cannot  be  aligned  and  a residual  cross-coupling 
will  result.  The  problem  is  completely  analogous  to  aligning  individual 
components  to  navigation  axes. 

Since  traditional  components  are  usually  misaligned  with  respect 
to  the  navigation  axes,  the  compensation  requirements  of  a system 
comprised  of  multisensors  is  not  significantly  different  from  that  of 
a system  built  of  other  components. 

Magnetic  interactions  such  as  those  occurring  in  the  feasibility 
multisensor  cannot  be  distinguished  from  a magnetic  misalignment. 

Because  the  gyro  torquer  is  mounted  on  the  magnetic  structure  shared 
with  the  radial  transducer,  case  rates  about  the  input  axis  cause  a 
shift  in  the  radial  null.  Because  of  elastic  restraint  a force  will 
be  read  on  the  accelerometer  channels.  Thus  the  relation  between  case 
rate  and  radial  torque  appears  as  a cross-coupling  coefficient. 

3.4  THERMALLY  INDUCED  ERRORS 

The  acceleration  measurement  of  the  multisensor  is  the  torque 
required  to  balance  the  couple  established  by  the  float's  centers  of 
gravity  and  buoyancy,  a fact  which  is  developed  in  Appendix  A and  in 
Appendix  F.  If  the  magnitude  or  position  of  the  equivalent  buoyant 
force  changes,  erroneous  acceleration  readings  ensue. 

As  detailed  in  Appendix  E,  a temperature  nrauient  along  the 
output  axis  can  cause  the  location  of  the  buoyant  force  to  change 
while  the  magnitude  remains  constant.  Because  fluid  motion  and  shear 
torques  accompany  this  effect,  the  term  convection  torque  is  also 
used.  The  analysis  of  Appendix  E predicts  a torque  sensitivity  of  4 
dyn-cm/g°F  (change  from  end  to  end  along  the  output  axis) . With  pen- 
dulosity  of  0.8  gm-cm,  the  scale  factor  sensitivity  is  5000  ppm/°F. 

If  the  temperature  level  changes  uniformly,  the  magnitude  of  the 
buoyancy  coefficient  will  change  because  of  the  flotation  fluid's 
sensitivity  to  temperature.  The  change  in  flotation  force  is  balanced 
by  the  radial  suspensions  which  are  aligned  with  the  gravity  vector. 

If  torques  are  summed  about  the  point  through  which  the  buoyant  force 
acts,  no  change  is  read  in  the  torque.  This  cancellation  is  equivalent 


to  selecting  the  weighting  between  current  readings  from  the  two  sus- 
pensions so  that 


k = q34/q44 


in  (F-14) . Because  the  buoyancy  is  typically  very  sensitive  to  tem- 
perature (400  ppm/°F  for  the  feasibility  instrument)  , removal  of  the 
buoyancy  level  change  from  the  torque  summation  is  desired. 

The  multisensor  feasibility  instrument  was  designed  symmetrically 
from  end  to  end  so  that  temperature  gradients  would  be  minimized.  A 
temperature  control  stabilized  level  changes  and  the  current  readings 
from  the  suspensions  at  either  end  of  an  axis  were  tuned  to  eliminate 
the  buoyancy  level  effects  on  the  torque. 

Assuming  temperature  gradient  controlled  to  0.01°F,  a scale  factor 
change  of  50-100  ppm  is  expected.  With  a float  mass  of  140  gm,  a fluid 
density  sensitivity  of  400  ppm/°F,  0.8  gm-cm  pendulosity,  5.6  cm  between 
suspensions,  temperature  level  change  of  0.01°F,  and  the  scalinq  between 
ends  correct  to  3 percent,  the  scale  factor  changes  57  ppm.  These  values 
are  approximately  consistent  with  those  seen  in  the  testing.  However, 
uncertainties  in  null  position  could  contribute  to  scale  factor  error 
as  discussed  in  paragraph  3.6.  Because  of  the  large  torque  capability 
required  for  300°/s  slewinq  and  because  of  the  excellent  gyro  perfor- 
mance obtained,  better  scale  factor  performance  can  be  obtained  by 
increasing  the  pendulosity. 


3.5  ACCELEROMETER  ERRORS  INTRODUCED  BY  WHEEL  ANISOELASTICITY  AND 

WHEEL  HUNT 

The  rotational  equations  of  motion  (3-4)  through  (3-6)  con- 
tain torques  which  are  sensitive  to  acceleration  products.  These 
terms  reflect  compliance  of  the  float  structure  and  the  wheel  bearings. 
The  acceleration  products  are  presently  compensated  in  single-degree- 
of-freedom  gyros,  in  pendulous  accelerometers  and  in  specific  force 
integrating  receivers;  thus,  these  terms  were  only  considered  in  the 
axes  location  tests,  with  proper  compensation,  the  acceleration 
products  should  have  only  small  effect  on  multisensor  performance. 

The  motion  about  the  spin  axis  depends  on  the  wheel's  hunting, 
the  H term  in  (3-4) . The  magnitude  of  H may  be  estimated  by  observing 
the  ripple  on  the  wheel  power.  Typical  numbers  for  muitisensor  feasi- 
bility size  instruments  with  hysteresis  motors  are  10  milliwatts 
ripple,  peak-to-peak  at  a frequency  of  10  Hz.  With  a wheel  speed  of 
48,000  r/min  (the  feasibility  instrument's  speed),  the  rms  torque  ripple 
is  7 dyn-cm  rms  at  i0  Hz.  With  800  gm-cm  pendulosity,  the  rms  velocity 
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error,  0.14  cm/s,  is  less  than  the  velocity  accuracy  of  0.7  cm/s 
required  for  a 50  ppm  system  as  estimated  in  Appendix  D. 

3.6  ACCELEROMETER  RESOLUTION  CONSIDERING  FORCER  ELASTIC  RESTRAINT 

AND  SIGNAL  GENERATOR  NOISE  AND  NULL  SHIFT 

Based  on  experimental  observations  and  the  analysis  discussed  in 
Appendix  D,  the  effects  of  the  accelerometer's  elastic  restraint  and 
the  accelerometer  signal  generator's  noise  upon  the  rebalance  force  is 
discussed. 

Elastic  restraint  is  defined  as  the  partial  deviation  of  the 
suspension's  force  with  respect  to  displacement.  An  ideal  radial 
suspension  forcer  would  exert  a force  which  depends  only  on  current; 
however,  normal  forcers  also  depend  on  radial  position.  Elastic 
restraint  must  be  considered  in  multisensor  design  for  the  following 
reasons : 

(1)  Destabilizing  influence  on  control  loops; 

(2)  Velocity  error  for  acceleration  step  change;  and 

(3)  Erroneous  force  reading  with  accelerometer  signal  generator 
null  shift. 

The  feasibility  instrument  forcer  primary  currents  are  fixed 
amplitude  square  waves  whose  measured  elastic  restraint  was  3 n 10' 
dyn-cm  per  end.  With  proper  tuning  of  the  signal  generator  as  a 
passive  suspension,  the  elastic  restraint  can  be  reduced  to  1 percent 
of  the  value  obtained  from  the  forcer  alone.  The  1 percent  figure  is 
conservatively  based  on  0.U01  cm  displacement  during  calibration  and 
ability  to  discern  torques  of  0.08  dyn-cm. 

For  typical  damping  constants  and  loop  gains  the  controller 
gains  are  sufficiently  greater  than  the  elastic  restraint  so  that  stability 
and  pole  placement  are  not  affected  even  if  the  SG  is  not  compensated. 

The  elastic  restraint  affects  numerator  dynamics  and  is  discussed 
under  velocity  error  and  signal  generator  null  shift. 

For  a Btep  change  in  acceleration,  the  steady-state  velocity  error 
is  proportional  to  the  effective  elastic  restraint.  With  compensation, 
the  steady  velocity  error  for  a one  gravity  step  change  is  0.006  cm/s 
compared  to  the  0.7  cm/s  commensurate  to  50  ppm  navigation  (see 
Appendix  D) . Thus,  with  compensation,  elastic  restraint  is  insignifi- 
cant and  without  compensation,  borderline. 
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If  the  signal  generator's  null  shifts,  the  float’s  will  be 
positioned  so  that  the  SG  output  is  zero.  With  elastic  restraint,  the 
position  change  will  cause  the  current  to  change  so  that  the  force  on 
the  float  remains  unchanged.  This  current  change  will  be  interpreted 
as  a bias  error.  If  the  null  shift  occurs  during  calibration,  it  can 
also  affect  the  scale  factor  determination. 

with  series  parallel  connection  used  in  the  feasibility  device's 
signal  generator,  the  null  is  directly  proportional  to  the  difference 
in  the  capacitance  of  the  parallel  legs.  As  calculated  in  Appendix  D, 
capacitance  stability  of  the  order  of  25  ppm  or  one  picofarad  is 
required  for  the  bias  stability  of  100  ug  in  the  feasibility  instrument 
which  was  not  accurately  tuned.  If  the  capacitors  shift  equally,  the  null 
and  the  force  reading  are  unaffected. 

The  capacitors  used  in  the  feasibility  device  had  a temperature 
sensitivity  of  approximately  200  ppm/0!'  and  the  wires  from  float  to 
capacitors  were  strapped  to  minimize  capacitance  changes . Despite  the 
multisensor ' s satisfying  the  bias  stability  requirement,  the  rotation 
during  the  tests  could  affect  the  capacitance  temperatures  and  appear 
in  the  scale  factor  stability  since  one  assumes  that  scale  factor  and 
bias  do  not  change  during  the  two  point  test. 

The  null's  sensitivity  to  capacitance  can  be  almost  eliminated 
by  connecting  the  signal  generator  windings  in  a bridge  with  capacitor 
across  the  bridge  for  tuning.  Since  null  is  determined  by  the  voltage 
across  the  capacitor,  capacitance  shifts  affect  null  position  only 
slightly  (depending  on  the  asymmetry  of  the  windings) . signal  genera- 
tor tuning  to  eliminate  elastic  restraint  and  better  suspension  capaci- 
tors (20  ppm/°F  can  be  obtained)  should  be  evaluated  before  additional 
end  housing  taps  required  for  bridging  are  inserted  into  the  multi- 
sensor  design. 

In  the  previous  single-degree-of-freedom  devices,  measurements 
were  taken  about  the  output  or  symmetric  axis  where  the  damping  con- 
stant is  typically  2 < 10s  dyn-cm-s.  Damping  constants  about  trans- 
verse axes  are  three  to  four  orders  of  magnitude  higher  because  of  film 
squeeze.  For  a given  bandwidth,  the  electronic  gain  is  proportional  to 
the  damping  as  detailed  in  Appendix  D;  thus,  small  amounts  of  pure 
signal  generator  noise  result  in  significant  rebalance  torque  errors. 

Pure  signal  generator  noise  is  defined  as  not  originating  in  actual 
motion  of  the  float.  Possible  sources  of  pure  siqnal  generator  noise 
are  Barkhausen  noise  in  the  magnetic  materials  and  preamplif ier  noise 
at  the  carrier  frequency. 
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If  the  signal  generator's  null  shifts,  the  float's  will  be 
positioned  so  that  the  SG  output  is  zero.  With  elastic  restraint,  the 
position  change  will  cause  the  current  to  change  so  that  the  force  on 
the  float  remains  unchanged.  This  current  change  will  be  interpreted 
as  a bias  error,  if  the  null  shift  occurs  during  calibration,  it  can 
also  affect  the  scale  factor  determination. 

with  series  parallel  connection  used  in  the  feasibility  device's 
signal  generator,  the  null  is  directly  proportional  to  the  difference 
in  the  capacitance  of  the  parallel  legs.  As  calculated  in  Appendix  r, 
capacitance  stability  of  the  order  of  25  ppm  or  one  picofarad  is 
required  for  the  bias  stability  of  100  ug  in  the  feasibility  instrument 
which  v/as  not  accurately  tuned.  If  the  capacitors  shift  equally,  the  null 
and  the  force  reading  are  unaffected. 

The  capacitors  used  in  the  feasibility  device  had  a temperature 
sensitivity  of  approximately  200  ppm/°F  and  the  wires  tr  float  to 
capacitors  were  strapped  to  minimize  capacitance  changes.  Despite  the 
multisensor's  satisfying  the  bias  stability  requirement,  the  rotation 
during  the  tests  could  affect  the  capacitance  temperatures  and  appear 
in  the  scale  factor  stability  since  one  assumes  that  scale  factor  and 
bias  do  not  change  during  the  two  point  test. 

The  null's  sensitivity  to  capacitance  can  be  almost  eliminated 
by  connecting  the  signal  generator  windings  in  a bridge  with  capacitor 
across  the  bridge  for  tuning.  Since  null  is  determined  by  the  voltage 
across  the  capacitor,  capacitance  shifts  affect  null  position  only 
slightly  (depending  on  the  asymmetry  of  the  windings) . Signal  genera- 
tor tuning  to  eliminate  elastic  restraint  and  better  suspension  capaci- 
tors (20  ppm/°F  can  be  obtained)  should  be  evaluated  before  additional 
end  housing  taps  required  for  bridging  are  inserted  into  the  multi- 
sensor design. 

In  the  previous  single-degree-of-freedom  devices,  measurements 
were  taken  about  the  output  or  symmetric  axis  where  the  damping  con- 
stant is  typically  2 * 10 5 dyn-cm-s.  Damping  constants  about  trans- 
verse axes  are  three  to  four  orders  of  magnitude  higher  because  of  film 
squeeze.  For  a given  bandwidth,  the  electronic  gain  is  proportional  to 
the  damping  as  detailed  in  Appendix  D;  thus,  small  amounts  of  pure 
signal  generator  noise  result  in  significant  rebalance  torque  errors. 

Pure  signal  generator  noise  is  defined  as  not  originating  in  actual 
motion  of  the  float.  Possible  sources  of  pure  signal  generator  noise 
are  Barkhausen  noise  in  the  magnetic  materials  and  preamplif ier  noise 
at  the  carrier  frequency. 
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An  alternate  viewpoint  argues  that  for  a given  loop  bandwidth 
increasing  damping  decreases  the  displacement  for  a given  input  so  that 
eventually  the  actual  motions  which  determine  the  forces  are  lost  in  the 
signal  generator  noise. 

For  50  ppm  navigation  which  implies  a velocity  error-of-error  of 
0.7  cm/sec  and  to  expedite  calibration  so  that  filtering  and  excessive 
sampling  (10,000)  points  are  not  required,  a damping  constant  less  than 
10 7 dyn-cm-s  is  required, 

3.7  REACTION  TORQUES  CAUSED  BY  AXIAL  SUSPENSION 

Previous  multisensor  testing  demonstrated  a large  reaction  torque 
in  the  transverse  channels  caused  by  the  axial  suspensions.  With 
passive  suspensions  a quiescent  force  is  exerted  at  each  end  which 
oppose  one  another.  With  machining  tolerances,  these  forces  are  not 
colinear  so  that  a couple  exists  about  a transverse  axis.  The  multi- 
sensor was  wired  with  active  time  shared  axial  suspensions  so  that  the 
quiescent  forces  with  perfect  flotation  would  be  zero,  and  hence,  the 
reaction  torques  would  be  reduced. 

3.8  ADJUSTMENTS  OF  DAMPING  CONSTANTS— INTERACTION  BETWEEN  INPUT  AND 

OUTPUT  AXES  COUPLED  THROUGH  ANGULAR  MOMENTUM 

The  multisensor's  three  dimensional  damping  characteristics 
should  be  altered  for  several  reasons.  As  discussed  in  Section  3.7, 
deunping  about  the  transverse  axes  must  be  reduced  to  decrease  noise 
in  the  torque  rebalance  signal  caused  by  signal  generator  uncertainty. 

A decrease  in  radial  and  axial  damping  constants  would  decrease  reac- 
tion time  because  of  more  rapid  transients.  Lower  fluid  viscosity 
facilitates  filling  and  reduces  contamination  and  density  gradient 
problems . 

In  reducing  fluid  viscosity  and  damping  constants,  the  following 
effects  must  be  considered! 

(1)  Time  constants  for  single  axes;  and 

(2)  Two-degrees-of-freedom  motion  coupled  through  wheel  angular 
momentum. 

For  this  discussion,  a single-degree-of-freedom  gyroscope  (or, 
alternatively,  a gyro  which  obeys  single-axis  theory)  may  be  controlled 
by  designing  the  torques  about  the  output  axis  as  a dynamic  function  of 
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the  float  rotation  about  the  output  axis.  Similarly,  the  input  axis 
control  torques  are  functions  of  the  float-case  angle  about  the  input 
axis. 

When  reducing  damping,  single  axis  effects  such  as  the  time  con- 
stants determined  by  the  damping  constant  divided  by  the  moment  of 
inertia  and  the  permissible  float  motions  must  be  considered.  Since 
the  design  proposed  here  will  result  in  time  constants  equivalent  to 
those  presently  used,  this  topic  will  not  be  pursued. 

When  altering  damping  constants,  the  interaction  of  input  and 
output  axis  loops  which  are  coupled  through  the  wheel  momentum  must  be 
considered.  With  active  suspensions  whose  bandwidths  are  comparable 
to  those  of  multisensor  applications,  the  single  axis  theory  is  valid 
(the  poles  of  the  coupled  system  approximate  those  developed  from 
single-axis  theory)  when  the  product  of  the  damping  constants  about  the 
input  (Cj)  and  output  (CQ)  axis  is  greater  than  100  times  the  angular 
momentum  (H)  squared.  For  lower  damping  products,  suitable  control 
requires  that  the  torque  about  the  output  and  input  axes  be  functions 
of  both  the  input  and  output  angles. 

For  many  practical  designs  and  bandwidths.  Appendix  C demonstrates 
that  inertia  terms  may  be  negligible  even  with  low  damping.  Thus,  the 
motion  equations  (C-3)  and  (C-4)  for  the  input  axis  are  virtually 
identical  with  that  of  the  output  axis.  The  form  of  the  two  equations 
is  identical  except  that  the  coefficients  of  the  float-ccse  angles  are 
defined  by  different  quantities. 

For  present  Draper  designs, 

Cj  CQ  > 100  H2 

/ 

and  the  torque  about  the  output  axis  is  controlled  as  a function  of  the 
output-axis  angle  and  similarly  for  the  input  axis.  The  similarity  of 
the  two  equations  then  dictates  that  for: 

H2  < 100  Cj  CQ 

The  torque  about  the  input  axiB  should  be  a function  of  the  angle  about 
the  output  axis  and  similarly  the  torque  about  the  output  axis  should 
be  a function  of  the  angle  about  the  input  axis.  As  in  the  single  axis 
situation,  the  two  equations  are  effectively  decoupled;  however,  the 
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dynamics  of  the  individual  loops  are  now  cominated  by  the  angular  momentum. 
This  situation  is  identical  to  that  of  single-axis  theory  where  the  loop 
response  is  dominated  by  the  damping  constant. 

The  above  observation  may  be  applied  directly  to  multisensor  design. 
Consider  the  feasibility  instrument,  fluid  viscosity  could  be  selected 
so  that  the  damping  constant  about  the  spin  axis  is  roughly  2 * 10! 
dyn-cm-s,  a figure  which  allows  good  loop  response  and  little  problem 
from  signal  generator  noise.  The  input  and  output  deunping  product 
would  be  less  than  H2/100.  The  crossfeedback  described  in  the  previous 
paragraph  would  stabilize  the  input  and  output  axis  loops  with  excellent 
dynamic  characteristics.  That  is,  the  three  rotational  loops  could  be 
designed  with  suitable  dynamic  characteristics. 

3.9  CONTROL  OF  ROTATIONAL  AND  TRANSLATIONAL  MOTIONS 

With  cylindrical  multisensor  designs,  the  radial  positions  at 
both  ends  of  the  float  are  measured  and  radial  forces  are  applied  at 
both  ends.  Since  the  radial  forcers  are  multisensor  outputs,  the 
implications  of  forcer  connections  and  control  design  are  discussed. 

The  radial  position  signals  indicate  the  displacement  of  the 
float's  ends  with  respect  to  the  null  position.  Because  the  permitted 
float-to-case  rotations  are  small,  the  signals  can  be  summed  to  yield 
translation  of  the  float  center  or  differenced  to  obtain  float  rota- 
tion. 

With  double  wound  forcers  which  are  included  in  the  feasibility 
instrument,  the  following  design  options  for  forcer  connections  are 
defined t 

(1)  End  control  (Figure  3-1)  where  the  forces  in  each  end 
are  controlled  independent ly»  and 

(2)  Mode  control  (Figure  3-2)  where  the  windings  are  connected 
so  that  ideally  forces  and  torques  about  the  float  center 
may  be  applied  independently. 

The  mathematical  implications  of  these  controls  are  discussed  in 
Appendix  B. 

Presently  the  feasibility  multisensor  exploys  end  control  with 
the  force  at  a given  end  controlled  by  the  displacement  at  that  end. 

This  option  was  selected  because  of  the  ease  with  which  the  precision 
pattern  field  could  be  obtained. 
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Figure  3-1.  Schematic  of  radial  forcer  connections, 
individual  end  control. 

The  dynamic  compensators  which  related  position  to  current 
(hence,  force)  were  selected  to  obtain  the  desired  rotational  response 
while  accepting  an  underdamped,  more  slowly  decaying,  translational 
signature  as  discussed  in  Appendix  B.  The  slower  mode  did  not  hamper 
the  feasibility  instrument's  performance;  however,  independent  loop 
tuning  can  be  obtained  by  crossfeedback,  were  the  displacement  at  the 
opposite  end  is  also  used  to  control  the  force.  Crossfeedback  can  be 
used  to  electronically  decouple  the  modes  as  discussed  in  Appendix  B. 

The  common  mode  rejection  of  the  translational  motion  (which  is 
undesirable  because  of  its  strong  temperature  dependence  as  discussed 
in  paragraph  3.4)  was  accomplished  in  calculations  external  to  the 
control  loops.  Because  of  the  uncertainty  in  establishing  the  scale 
factors  of  the  two  forcers  (Appendix  F and  Section  2)  , probably  only 
95  - 97  percent  of  the  common  mode  was  rejected. 
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Figure  3-2.  Schematic  of  radial  forcer  connection, 
separate  rotational  and  translational 
mode  control. 

The  mode  control  permits  the  common  mode  to  be  rejected  mechani- 
cally. Because  of  forcer  and  instrument  dissimilarity,  the  rejection 
would  probably  be  97  percent  successful.  But  the  rotational  torque 
could  be  trimmed  by  adding  a small  portion  of  the  translational  force 
so  that  mode  control  promises  common-mode  rejection  greater  than  99 
percent. 

The  separation  of  the  modes  is  required  to  implement  the  lightly 
damped  mechanization  described  in  the  previous  section.  Once  the 
rotational  torques  are  separated  from  the  forces,  the  connections  can 
be  arranged  to  permit  single-degree  or  lightly  damped  operation. 

3.10  MAXIMUM  CASE  RATES 

The  work  statement  specifies  that  case  rates  of  300  deg/s  (or 
5 rad/s)  should  be  considered.  In  determining  maximum  permissible 
case  rate,  the  following  factors  must  be  considered! 
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(1)  At  steady-state  is  sufficient  torque  available  about  both 
the  input  and  output  axes  to  rebalance  the  momentum-case 
rate  product; 

(2)  Provided  sufficient  torque  is  available,  can  the  damping  and 
loop  bandwidth  be  tuned  so  that  the  float  does  not  hit  the 
stops  during  a transient;  and 

(3)  Will  terms  such  as  the  inertia  difference  times  the  product 
of  case  rates  require  compensation  in  navigation. 

Appendices  B,  C,  or  D may  be  extrapolated  to  demonstrate  that  the 
float  can  be  restrained  from  case  contact  with  the  proposed  integral 
control  loops.  Assume  a change  from  0 to  5 rad/s  can  be  represented  by 
a ramp  of  1 rad/s/s.  From  (D-4) , for  the  22  Hz  bandwidth,  tuning  with 
angular  momentum  equal  loop  damping,  the  displacement  is  600  pin. 
However,  the  problem  of  resolving  meru  with  maximum  rates  of  1 rad/s 
must  be  resolved  (see  Section  4) . 

For  strapdown  application,  compensation  for  the  cross-product 
terms  will  be  required.  An  angular  momentum  of  20,000  dyn-cm  corres- 
ponds to  0.0015  dyn-cm.  The  difference  between  moments  of  inertia 
about  spin  (approximately  equal  input)  and  output  axis  is  roughly 
5 gm-cm  for  a 10-size  device.  At  a slew  of  0.1  rad/s,  the  inertia 
difference  results  in  a torque  of  0.05  dyn-cm  so  that  compensation  is 
required. 
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SECTION  4 


DESIGN  PROJECTION 

A preliminary  design  of  a 10-size  multisensor  (nominal  case  diam- 
eter of  one  inch)  appears  in  Figure  4-1.  Important  design  parameters 
are  listed  in  Table  4-1.  A brief  discussion  of  the  design  highlights 
is  included. 

The  10-size  multisensor  is  designed  with  pendulosity  of  1 gm-cm, 
slightly  larger  than  the  feasibility  instruments.  At  10  gravity,  the 
required  rebalance  torque  is  10,000  dyn-cm,  a figure  within  the  forcer 
capability  by  a factor  of  four.  Since  thermally  induced  errors  scale 
as  a characteristic  dimension  to  the  fourth  power  (paragraph  31.4,  and 
Appendix  E) , a factor  of  ten  improvement  in  acceleration  sensitive 
stability  could  be  realized. 

The  10-size  PIG  (Pendulous  Integrating  Gyroscope)  has  demonstrated 
torque  uncertainties  which  would  give  performance  better  than  the  one 
meru  design  guidelines  with  either  the  5,000  or  10,000  dyn-cm  wheel. 

To  reduce  compliance  torques,  hemispherical  gas  bearing  wheels  are 
selected.  Although  hysteresis  motors  perform  acceptably  (paragraph  3.5), 
further  reduction  of  wheel  hunting  and  more  stable  thermal  environment 
are  products  of  the  permanent  magnet  wheel  motor  design  depicted  in 
Figure  4-1. 

While  smaller  values  of  angular  momentum  could  satisfy  perfor- 
mance, 5000  were  selected  so  that  the  gyro  axes  would  operate  in  the 
lightly  damped  mode  where 


CI  °0  < 


H2 

100 


In  the  lightly  damped  mode  all  axes  could  be  designed  with 
similar  dynamic  response  as  described  in  paragraph  3.3.  The  rotational 
and  translational  pull-in  times  are  less  than  9 seconds  so  that  excel- 
lent warm-up  characteristics  will  be  achieved.  The  low  viscosity  fluid 
will  greatly  facilitate  fill  procedure. 
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Figure  4-1.  Multisensor  (low  damping  - 10  size). 


TABLE  4-1.  DESIGN  PARAMETERS  FOR  10-SIZE  MULTISENSOR. 


Float  Mass 

22  gms 

Moments  of  Inertia  About  S 

20  gm-cm2 

Moments  of  Inertia  About  I 

20  gm-cm2 

Moments  of  Inertia  About  0 

15  gm-cm2 

Angular  Momentum 

5 x 105  — 2 x 101*  dyn-cm 

Pendulosity 

1 gm-cm 

Fluid  Viscosity  w 

3 cP 

Deunping  Constants 

Radial  Rotation  (CR) 

3 x 10"  cP 

OA  Rotation  (CQ) 

30  cP 

Radial  Translation  (C,j,) 

1.5  x 10 5 cP 

Axial  Translation  (C^) 

1.5  x 10s  cP 

Torque  Maximum 

Output  Axis 

5 - 10,000  dyn-cm 

Transverse 

4 x 10**  dyn-cm 

Axial  motion  and  the  inherent  coupling  with  the  rotational 
measurement  axis  is  reduced  by  the  outward  facinq  mass  compensated 
bellows  which  the  feasibility  instrument  did  not  possess. 

To  mechanize  the  control  required  for  light  damping,  the  forcer 
windings  are  double  wound  and  connected  in  the  mode  control  of  paragraph 
3.9.  As  demonstrated  in  the  feasibility  instrument,  active  control 
will  be  used  to  adjust  dynamic  control  and  to  assure  that  the  float's 
position  with  respect  to  the  case  remains  constant.  The  radial  and 
axial  suspension  rotors  overlap  their  respective  stators  so  that 
cross-coupling  between  axes  is  minimized.  Because  of  the  lightly- 
viscous  damping  fluid,  the  magnetic  suspension  clearances  can  be  set 
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to  optimize  suspension  force  and  efficiency.  With  high  deunping,  reduc- 
tion of  the  transverse-axis  damping  resulted  in  increased  suspension 
clearance  and,  hence,  reduced  magnetic  efficiency. 

With  permanent  magnet  gyro  torquer  with  5000  dyn-cm  capability, 
case  rates  of  1 - 2 rad/s  could  be  balanced.  Since  current  readouts 
can  only  scan  six  orders  of  magnitude  and  since  the  ratio  of  1 rad/s  to 
1 meru  is  1.4  x 107,  a design  trade-off  between  maximum  case  rate  and 
resolution  must  be  made.  The  following  options  are  being  considered 
for  generating  the  radial  signals  required  for  control: 

(1)  Two  permanent  magnet  torquers  about  the  output  axis  with 
capacitive  pick-offs  in  the  annulus  to  generate  both  radial 
and  output  axis  signals.  This  option  is  shown  in  Figure  4-1. 
Because  of  the  instrument's  small  size,  questions  over  gyro 
signal  resolution  must  be  resolved.  Because  of  their  high 
impedance  and  the  goal  of  not  requiring  additional  flex 
leads,  the  capacitor  signal  generators  require  careful 
design  to  maintain  elastic  restraint  at  acceptable  levels 
for  gyro  signal  generator. 

(2)  One  permanent  magnet,  a microsyn  signal  generator,  and 
capacitor  pick-offs  for  the  radial  signals.  The  gyro  would 
be  controlled  as  in  present  instruments  where  the  microsyn 
signal  generator  has  proven  itself.  The  capacitors  could 
then  be  used  in  their  most  sensitive  normal  operation  to 
generate  the  radial  signals. 

(3)  Use  soft  iron  torquer  and  microsyn  which  double  as  the 
radial  signal  generator  as  in  the  feasibility  instrument. 

With  the  low  viscosity  fluid,  the  noise  propagation  problem 
which  appeared  in  the  feasibility  instrument  can  be  avoided. 
The  cross-coupling  between  forcer  axes  and  gyro  torquer  can 
be  compensated  (the  required  term  will  exist  in  any  case) 

as  described  in  paragraph  3.3.  Although  the  maximum  torques 
generated  with  soft  iron  are  an  order  of  magnitude  less  than 
those  of  the  permanent  magnet  torquer,  present  read-out 
devices  can  not  scan  the  seven  orders  of  magnitude  between 
1 rad/s  and  1 meru. 

(4)  Two  microsyns  and  magnetic  torquer.  The  microsyns  could  be 
used  to  simultaneously  measure  angle  about  the  output  axis 
and  radial  position  as  done  in  the  present  multisensor.  To 
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minimize  cross-coupling  a separate  torquer,  either  permanent 
magnet  or  soft  iron,  is  added,  a torquer  or  signal  genera- 
tor could  be  added  about  the  instrument's  circumference 
(belly)  or  added  in  the  end  housing  areas.  These  options 
require  further  refinement. 


SECTION  5 


CONCLUSIONS  AND  RECOMMENDATIONS 


5.1  CONCLUSIONS 

Based  on  test  results  of  a feasibility  instrument,  analyses,  and 
design  studies,  the  multisensor  concept  is  suitable  for  navigation  and 
should  be  pursued  further. 

With  an  instrument  which  was  neither  designed  nor  optimized  as  a 
multisensor,  accelerometer  stability  of  50  \ig  bias  and  75  ppm  scale 
factor  were  achieved  in  two  position  100-hour  tests.  Despite  the 
large  (by  gyro  standards)  output  axis  pendulosity,  the  gyro  exhibited 
bias  stabilty  of  0.1  meru  and  acceleration  sensitive  stability  of  0.3 
meru/g. 

Because  of  high  damping  about  the  transverse  axes,  and  uncom- 
pensated bellows,  transient  behavior  was  not  at  operational  standards; 
however,  the  sources  of  the  transients  are  understood  and  design  altera- 
tions have  been  presented  which  circumvent  the  deficiencies  (see 
Section  4) . The  projected  size  10  design  will  have  radial  and  axial 
pull-in  times  less  than  one  second  because  of  the  low  viscosity  fluid 
which  is  possible  because  of  the  cross-connected  input-output  axis 
control.  By  increasing  pendulosity  and  because  of  the  temperature 
dependent  error's  sensitivity  to  size,  the  projected  design  should  have 
accelerometer  performance  significantly  better  than  the  feasibility 
instrument. 

The  multisensor  program  is  more  significant  than  a mere  demon- 
stration of  a multisensor  which  can  simultaneously  measure  angular  rate 
and  acceleration.  Most  importantly,  torque  summation  about  the  trans- 
verse axis  of  a cylinder  has  been  demonstrated.  With  the  ability  to 
sum  torques  about  a transverse  axis,  one  can  envision  an  entire  family 
of  inexpensive  multifunction  instruments.  From  identical  hardware  with 
only  slight  modifications  to  the  float,  multisensors,  two-degree-of- 
freedom  gyros,  two-degree-of- freedom  accelerometers,  and  three-degree- 
of -freedom  angular  rate  sensors  can  be  built. 
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5.2  RECOMMENDATIONS 

Recommendations  for  further  testing  of  the  present  instrument  and 
for  design  activities  are  listed. 

To  demonstrate  improved  performance,  the  present  feasibility 
instrument  should  be  tested  with  the  following  changes: 

(1)  More  insight  into  accelerometer  bias  could  be  obtained  by 
replacing  the  present  200  ppm  suspension  capacitors  with  20 
ppm  capacitors. 

(2)  Accelerometer  bias  stability  improvements  from  more  accurate 
cancellation  of  the  accelerometer  elastic  restraint  should 
be  investigated. 

The  feasibility  multisensor  should  be  tested  in  continuous  tumble 
tests.  The  tumble  test  will  allow  the  multisensor  to  be  evaluated  in  a 
dynamic  environment  and  will  investigate  the  stability  of  the  axes. 

The  feasibility  instrument  should  be  rewired  externally  to  mechan- 
ize the  mode  control  of  paragraph  3.9.  The  output  axis  signal  generator 
will  control  the  output  axis  torques  and  similarly  for  the  input  axis. 
Stability  and  sensitivity  tests  would  then  be  conducted  in  order  to 
compare  performance  with  mode  control  to  that  with  end  control. 

If  mode  control  is  successfully  demonstrated,  the  instrument  should 
be  disassembled,  cleaned,  and  refilled  with  fluid  whose  viscosity  is 
3 cP.  The  control  would  be  rewired  so  that  the  output-axis  torque  is 
controlled  by  the  float-case  angle  about  the  input  axis.  The  lightly 
damped  tests  will  verify  the  instrument  design  projected  in  Section  4. 
Particular  points  to  be  demonstrated  are: 

(1)  Fast  transients  with  axial  and  radial  pull-in  times  of  order, 
one  second; 

(2)  Stability  with  lightly  viscous  damping  fluid; 

(3)  Similar  dynamic  performance  about  the  three  rotational 
axes; 

(4)  Reduction  of  noise  in  accelerometer  channels. 

In  conjunction  with  the  testing  described  above,  design  activity 
is  required  with  particular  emphasis  on  signal  generator  design  and 
operational  requirements.  As  discussed  in  Section  4,  at  least  four 
options  are  available  for  generating  the  output  axis  angle  and  radial 
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displacements  necessary  to  implement  active  control.  The  options  and 
trade-offs  should  be  carefully  studied  to  establish  a design  optimized 
with  respect  to  performance  and  cost. 

In  order  to  design  suitable  instruments,  the  multisensor's  system 
applications  should  be  studied  so  that  the  feasibility  programs  are 
geared  toward  the  proper  operating  environment.  In  particular,  maximum 
slew  rates  and  minimum  sensitivity  (paragraph  3.10  and  Section  4)  should 
be  resolved. 

While  the  projected  design  of  Section  4 was  based  on  beryllium 
technology,  design  effort  should  be  focused  on  low-cost  multisensors 
made  from  molded  components . 
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APPENDIX  A 


MULTISENSOR  EQUATIONS  OF  MOTION 


A. 1 SUMMARY 

Equations  of  motion  are  derived  for  the  multisensor.  The  float 
is  modeled  as  a rigid  body  which  has  six-degree-of-freedom  with  respect 
to  the  case.  The  wheel  is  modeled  as  a rigid  cylinder  which  spins  with 
respect  to  the  float  about  its  (the  wheel's)  symmetric  axis.  While 
several  simplifications  and  assumptions  are  made,  the  characteristics 
essential  for  multisensor  analysis  are  retained.  Because  several  of 
the  six-degrees-of-freedom  are  only  loosely  coupled,  the  six  equations 
of  motion  are  not  solved  simultaneously.  Additional  appendices  inves- 
tigate in  more  detail  certain  characteristics  of  the  equations  which 
directly  impact  multisensor  performance 

A. 2 NOMENCLATURE 

For  the  ensuing  analysis  the  following  conventions  will  be 
followed: 

(1)  Vectors  will  be  denoted  by  an  overline  and  matrices  by 
double  overline.  A superscript  indicates  the  reference 
frame  in  which  the  vector  is  measured.  Two  letter  sub- 
scripts will  be  employed.  The  fiist  subscript  indicates 
the  coordinate  frame  where  the  vector  originates  and  the 
second,  the  ending  point  of  the  vector.  For  example,  a^c 
indicates  the  acceleration  of  the  case  with  respect  to 
inertial  space  in  float  coordinates.  The  vector  triad  is 
ordered  as  spin,  input  and  output. 

(2)  The  acceleration  (a°c)  of  the  case's  nominal  center  and 
the  angular  rotation  (u>^c)  of  the  case  with  respect  to 
inertial  space  will  be  considered  as  inputs.  a^c  includes 
gravitational  attraction. 


(A-la) 


(A-lb) 


where  S,  I,  and  O denote  spin,  input,  and  output,  respec- 
tively. 

(3)  Three-degrees-of- freedom  are  given  by  displacement  vector 

from  the  nominal  center  of  the  case  to  the  nominal  center  of 
the  float;  that  is, 


(A-2) 


(4)  Three  additional  degrees  of  freedom  are  defined  by  the 
angular  rotation  of  the  float  with  respect  to  the  case. 

The  rotation  matrix  from  float  to  case  coordinates  is  given 
by; 


where  small  angles  have  been  assumed.  The  angles  are  measured 
about  the  nominal  case  axes  to  the  nominal  float  axes  as 
indicated  by  the  subscripts. 
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(5) 


The  acceleration  of  the  nominal  float  center  with  respect 
to  inertial  space  is  given  approximately  by: 


(A. 4a) 


where  the  products  of  case  rates  (ui?c)  and  float  to  case 
displacements  (r^p)  have  been  assumed  small. 

(6)  The  angular  velocity  of  the  float  with  respect  to  inertial 
space  in  float  coordinates  is  given  within  the  limit  of 
small  angle  approximations  by: 


for  this  use,  5^  can  be  assumed  identity, 
c 


(A. 4b) 


A.  3 EQUATIONS  OF  MOTION 

Application  of  Newton's  Laws  to  the  float  result  in  six  scalar 
euuations  which  permit  solution  for  the  six-degrees-of-freedom.  For 
linear  translation: 


and  for  rotation 


IF 


-cm 


(A- 5) 


(A-6) 


where 

E denotes  summation 
F - forces 

M ■ applied  torques 

H - angular  moment  of  the  float 
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- differentiation  with  respect  to  time 

m » mass  of  float  including  wheel 

a?  __  • acceleration  of  float  center  of  mass  which  is 

i-cm  _ 

approximately  a^p  by  design. 

Differentiation  of  the  float's  angular  momentum  is: 


dH 

3t 


d 

at 


[ ' =T,]  * “It  » [*!?,] 


d 

at 


Hw  + “ 


iF 


x H 


w 


(A-7) 


where : 


0 


0 


0 Ij  0 

O 0 i0 


is  the  inertia  matrix  or  sensor. 


I includes  both  the  inertias  of  the  float  and  the  wheel.  As  a first 
approximation,  cross  products  of  inertia  are  neglected. 


H. 


H 

0 

0 


spin  angular  momentum  of  wheel 


where : 

H " V 

Iw  ■ inertia  of  wheel  about  its  spin  axis 

fl  - angular  velocity  of  wheel  with  respect  to  float 


The  sum  of  the  forces  is  represented: 


EF 


Psus  + ^vis  + *buoy  + ^add 


(A-8) 
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where: 


Z.M  = 
where 


F ■ suspension  forces  which  are  discussed  further  in 
8 Appendices  B,  C,  and  F. 

F - fluid  forces  exerted  by  fluid  during  float  motion 


the  fluid  forces  are  assumed  to  functions  of  the  trans- 
lation velocity  in  the  given  direction;  however,  this 
simple  model  may  be  inaccurate.  Reference  5 discusses 
more  accurate  models. 

F.  = buoyant  forces  exerted  by  fluid.  The  buoyant  force  is 
” described  by  the  vector  relation: 


where  the  matrix  B is  a function  of  temperature  distri- 
bution. As  a first  approximation,  let  the  matrix  B equal 
a scalar  B multiplied  by  the  identity  matrix. 

The  sum  of  the  applied  torques  is  represented  as: 


sus 


+ M. 


vis 


+ M, 


buoy 


+ M 


V 


+ r_,  * F.  + M , , 

Fb  buoy  add 


(A-9) 


M * torques  applied  by  the  magnetic  suspensions  which  are 
8US  discussed  in  Appendices  B,  C,  and  F. 


Mv^a  * torques  applied  by  fluid  during  float  motion. 


As  for  the  viscous  forces,  this  simple  model  may  be 
inaccurate. 


105 


M 
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temperature  sensitive  torques  exerted  by  fluid  indepen- 
dent of  float  motion.  These  torques,  also  called  con- 
vection torques,  are  discussed  in  greater  detail  in 
Appendix  E. 

torques  caused  by  the  noncoincidence  of  the  center  of 
mass  and  the  origin  about  which  torques  are  summed.  Note 
that  the  nominal  center  for  torque  summation  is  arbitra- 
rily selected.  In  float  coordinates: 


rF-cm 


X 


where: 

_p 

rF-cm  ttie  displacement  from  the  origin  to  the  center  of 

mass,  a constant  vector  in  float  coordinates.  For  most 
multisensors,  study  a^F  can  be  set  equal  to  a^c  since 
the  float  case  angles  are  restricted  to  values  less  than 
0.001  rads  by  the  float  case  clearance.  The  input  of 
this  term  is  discussed  following  (A-ll) . 

M 2 - torques  proportional  to  acceleration  squared  which  may 

9 be  represented  as: 


"kiiaiao"kisasao“kioao2+koiai2+kosasai+kooaoai 


Mg*  “ | "kOOaOas"kOlalaS_kOSaS2+kSOa02+kSiaiaO+kSSaSaO 

'kssasarksoaoarksiai2+kisas2+kioaoas+knasai 


where: 

k.  . are  constant  coefficients  which  represent  structural 

13  and  bearing  compliance. 

r_.  x F.  is  caused  by  the  buoyant  force ' s not  acting  through 
fd  d oy  the  coordinate  center. 

M ..  are  other  torques  such  as  those  caused  by  flex  leads  or 
aaa  unmodeled  interactions  of  magnetic  fields. 

With  the  simplifying  assumptions  notes,  (A-7) , (A-3) , and  (A-9) 
are  substituted  into  (A-5)  and  (A-G)  to  obtain  the  equations  of  motion 
in  a form  useful  for  additional  study.  The  equations  of  motion  for 
translation  are: 
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S TS 


mX_  + CmXT 
I T I 


■*8  + CAXX0 


(B-ml  .s  ♦ f,u8 

(S)  + Fadd  (S) 

(A-lOa) 

<B-">  *1  * p.us 

(I)  + Fadd  (I) 

(A-lOb) 

IB-"°  fl0  + F,us 

(0)  + Fadd  (0) 

(A-lOc) 

For  rotation,  the  motion  equations  are: 

*S  0s  + CR  0S  = alP  a0  UI  + MSUS  (S)  + “g2 (S)  + Madd  (S) 

+ Mbuoy  (S)  " » " XS  “S  + (IfV(VV(VV 


(A-lla) 


XX  9I  + CR  6I  + H0O 


"aSP  + a0  US  + Msus  (I)  + Mg2(I)  + Madd  (I) 


+ Mbuoy  (I)  “ H“o  ‘ XI  “i  + <IQ-IS)  (W  (VV  (A'llb) 


X0  90  + c0  0o  ~ H9I  " as  aI  Us  + Msus  (0)  + Mg2  (O)  + Madd  (O) 


+ “buoy  (O)  + Ha)I  ' X0  “0  + (IS"II)  (W  (u,I +0I)  (A-llc) 


where: 

P ” mrF-cm(0)  ‘ BrFb  (0) 

- pendulosity  along  the  output  axis.  The  pendulous  torques 
which  the  suspensions  rebalance  are  the  couple  between  the 
acceleration  and  buoyant  forces. 

US  " mrF-cm(S)  " BrFb(S) 

■ unbalance  along  spin  axis. 

UI  ■ mrF-cm(I>  - Brpb(I) 

■ unbalance  along  input  axis. 

Ideally,  the  unbalance  terms  Ug  and  UI  should  be  zero. 
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APPENDIX  B 


COUPLED  ROTATIONAL  AND  TRANSLATIONAL  MOTION 


B.l  SUMMARY 

To  further  understanding  of  suspension  tuning  in  the  multisensor, 
a float  permitted  two-degree-of-freedom  (2-DOF)  is  studied.  The  per- 
mitted motion  is  rotation  and  linear  translation  with  forces  applied 
at  the  float  ends. 

The  control  requirements  necessary  to  decouple  the  modes  and  the 
effects  of  control  strategy  on  the  two  modes,  are  discussed.  The 
active  suspensions  used  for  the  feasibility  multisensor  are  studied. 

The  principal  conclusions  are: 

(1)  Float  suspensions  should  be  designed  with  consideration  of 
both  rotational  and  translational  modes. 

(2)  For  the  control  used  in  the  feasibility  instrument,  the 
rotational  mode  is  quicker  than  the  translational  mode. 

B.2  BASIC  MODEL 

As  depicted  in  Figure  B-l,  rotation  about  the  center  of 
mass  and  translational  motion  of  the  cm,  are  permitted.  The  equations 
of  motion  to  be  analyzed  are  Newton's  law  for  translational  motion: 

F + Fi  + F2  * x(ms2  + Cts)  (B-l) 

and  for  rotational  motion: 

M + F1t1  - F2*2  = 6 (Is2  + crs)  (B-2) 


where: 

s ■ the  Laplace  operator  d/dt. 

Fjy  f2  * suspension  forces. 

F - external  force  such  as  buoyancy. 
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B.3  EQUATIONS  OF  MOTION  WITH  SUSPENSION  FORCES 

Several  expressions  are  derived  to  relate  various  suspension  con- 
trol strategies  to  the  resulting  rotational  and  translational  inodes. 
General  expressions  are  derived  and  extraneous  terms  will  be  discarded 
when  appropriate. 

The  suspension  control  laws  written  in  terms  of  the  cm  motion  and 
rotation  about  the  cm  are: 


(B-3) 


where  the  G's  are  arbitrary  compensators  so  that  both  passive  and 
active  suspensions  can  be  handled.  In  matrix  form,  the  dynamic  equa- 
tions are: 


| G10+G20 

X 

*r 

F 

Glx*rG2x£2 

| IS’*cRs*GietrG2et2 

0 

M 

Since  signal  generators  usually  measure  the  displacement  at  the  float 
ends,  (B-3)  and  (B-4)  are  rewritten  with  x^  and  the  displacement, 
at  the  float  ends  (see  Figure  B-l) 


’ Xx  “I  [ 1 j If  X ' 

. *2  J " L 1_i  Jl  6 . 


(B-5) 


where: 

■ the  distance  from  end  1 to  the  cm  and  similarly  for  &2. 

M ■ externally  applied  torque  (for  example,  torques  caused  by 
fluid  motion) . 

m ■ float  mass. 

c^.  ■ damping  constant  for  translational  motion, 

x - displacement  of  center  of  mass. 
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I * moment  of  inertia  about  center  of  mass 
c ■ damping  constant  for  rotational  motion 
6 * angular  displacement  about  cm 

The  generalized  coordinates  x and  0 are  selected  since  these  vari- 
ables permit  the  simplest  description  of  the  resulting  motions,  speci- 
fically, if  pure  torques  are  applied  to  the  float,  the  resulting  motion 
can  be  viewed  as  a rotation  about  the  center  of  mass  which  does  not 
move  translationally . 

Solving  (B-5)  for  the  position  of  the  mass  center  and  the  angular 
position: 


(B— 6 ) 


Substitution  of  (B-G)  into  (B— 4 ) gives  the  suspension  control 
law  in  terms  of  the  end  displacements: 


(B— 7 ) 


This  work  will  discuss  strategies  for  controlling  the  suspension 
forces  in  general.  Specific  examples  pertaining  to  the  feasibility 
multisensor  will  be  presented. 

B.  4 DISCUSSION  ON  DECOUPLING 

With  freedom  to  select  the  four  compensators  of  (3)  , the  designer 
has  sufficient  degrees  of  freedom  to  decouple  the  translational  and 
rotational  modes  by  selecting: 

°26  ' -°ie 


Gix  l2  * Gie 


Gix  *1 


- G 


10 


1 


G2x  *2  + G20  ! 


G2x  *1 


- G 


20 


*1  + *2 
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(B-8b) 


G 


2x 


G 


lx 


Sufficient  freedom  then  remains  to  independently  tune  the  rotational 
and  translational  loop;  that  is,  with  (B-8)  , (B-4)  becomes: 


With  (B-8) , the  control  law  ..or  decoupled  motion  when  the  end  positions 
are  measured  is  determined  from  (B-7) : 


(B— 10) 


Thus,  in  general,  decoupled  motion  cannot  be  attained  by  controlling 
the  force  with  only  the  measurement  from  that  end.  Since  the  forces 
must  be  measured  at  both  ends,  the  complexity  of  mechanizing  (B-7)  is 
not  exorbitant.  In  summary,  the  requirements  for  complexity  decoupling 
the  loops  have  been  derived  in  (B-8)  and  (B-10) . 


B.  5 END  DISPLACEMENT  MEASUREMENT  WITHOUT  CROSSFEEDBACK 


The  multisensor  feasibility  is  tuned  with  active  suspensions  and 
compensators  G (s)  so  that: 


r 

r-l 

fc. 

1 

G (8) 

! 0 
i 

V 

L ?2 

m 

1 

1 o 

1 

1 

! G (S) 

X2 

L * J 

L 

1 

(B-ll) 
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A general  expression  for  systems  with  end  displacement  measurement  and 
without  crossfeedback  (F^  is  not  a function  of  x2  and  similarly  for  F2, 
Xj)  is  derived  from  (B-7)  by  setting: 


G 


10 


G . t. 
lx  1 


(B-12a) 


G 


20 


“G2x*2 


(B-12b) 


Substitution  of  (B-12)  into  (B-4)  yields  the  equivalent  of  motion  for 
the  rotational  and  translational  modes: 


' n,s’4CTS*Glx.G2x  | 
1 

Glxlll~G2xt2 

X 

F 

Glx£l"G2x*2  ! 

I,  +Cr.*GlxV«2xl2’ 

e 

M 

In  general,  the  modes  are  not  decoupled;  however,  for  most  practical 
systems  (even  the  multisensor  with  pendulosity  along  the  output  axis) 
is  approximately  and  Glx  is  approximately  G2x  so  that  the  trans- 
lational and  rotational  modes  are  effectively  decoupled.  Mass  m and 
inertia  I can  be  neglected  compared  to  the  damping  constant  CR  and  CT 
and  the  motion  equations  for  systems  which  measure  end  displacement 
and  do  not  employ  crossfeedback  become: 


■ CTs  + 2Glx  | 0 

1 

X " 

F " 

1 

1 

- 

0 i V + 2Glx*I  . 

0 

M 

Although  the  equations  are  decoupled,  the  freedom  to  independently 
tune  the  rotational  and  translational  loops  has  been  lost  (if  the 
lengths  are  assumed  fixed)  . This  loss  of  freedom  is  caused  by  the 
decision  not  to  employ  crossfeedback.  Thus,  if  one  loop  is  well  tuned, 
less  desirable  dynamic  performance  might  result  for  the  other  loops. 

With  passive  radial  suspensions,  both  modes  are  first  order  with 
different  characteristic  times.  For  active  suspension,  the  mode 
shapes  become  more  complicated  so  that  the  Glx<*>  selection  is  more 
interesting. 
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B. 6 THE  MULTISENSOR  RADIAL  SUSPENSION 


As  presently  configured,  the  multisensor  has  active  radial  sus- 
pensions that  result  in  dynamic  equations  of  the  form  (B-14) . This 
active  suspension  is: 


2G 


lx  1 


c„kR  (St.) 
S (S+b) 


(B-15) 


The  equations  of  motion  for  the  translational  and  rotational  modes  are 
then: 


where i 


e 


x 


[ 

[ 


S’+bS+kR(S+a) 


S,+bS+kT(S-»a) 


] 

] 


S(S+b) 

t'R 


(B-ICa) 


(B-16b) 


C_kB 
R R 

kT  cprp- 


Thus,  the  dynamic  performances  differ  because  of  the  different  loop 
gains  which  are  caused  by  the  rotational  and  translational  damping 
constants  not  being  equal  and  the  relation  between  force  and  torque 
through  the  moment  arms.  The  different  performance  are  then  related 
to  the  pole  positions  on  the  root  locus  diagram.  Selection  of  perfor- 
mance on  one  loop  decides  the  gain  and,  hence,  pole  positions  for  the 
other  loop. 

For  the  feasibility  multisensor: 


Cr/Ct  “ j (cm) 2 
- 2.8  cm 
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the  dynamic  compensator  was  originally  selected  so  that: 


1 

k_  - 1.41  x 10"  ~ 

R s' 

A = 22.8  s“ 1 

B = 173  s"1 


The  gain  of  the  translational  loop  is  then: 

kT  = 900  1/s2 


The  eigenvalues  for  the  rotational  mode  are: 
-34.4,  -69.7  + 67. 3j 


and  for  the  translational  mode: 

-168.,  -2.04  ± 10.2  j 

Further  insight  into  these  responses  is  offered  by  the  root  locus  of 
Figure  B-2.  Note  that  the  translational  loop  is  underdamped  (damping 
ratio  - 0.2)  and  decays  much  more  slowly  that  the  rotational  loop.  The 
responses  of  the  two  loops  are  compared  in  Figure  B-3  for  step  inputs  of 
force  and  torque. 

By  using  crossfeedback  (B-10) , the  two  loops  can  be  tuned  indepen- 
dently if  improved  performance  is  desired.  Note  that  the  loops  can  be 
decoupled  by  the  windings  discussed  in  paragraph  3.9.  The  cross-connection 
allows  application  of  the  translational  force  independently  of  the  rota- 
tional torque. 

Because  of  noise  described  in  Section  2 and  paragraph  3.6,  the  loops 
used  through  multisensor  tests  were  tuned  so  that  closed-loop  poles  were 
1/30  of  those  used  in  examples. 


115 


Figure  B.l.  Coordinates  for  rotational-translational  model 


APPENDIX  C 


TWO-DEGREE-OF-FREEDOM  FLOAT 
DYNAMICS  WITH  WHEEL  COUPLING 


C.l  SUMMARY 

The  behavior  of  a simplified  multisensor  model  which  focuses  on 
the  interaction  of  the  control  loops  about  input  and  output  axes 
coupled  by  the  wheel  angular  momentum  is  studied.  The  study  indicates 
the  conditions  under  which  control  loops  can  be  designed  from  single- 
loop analyses  and  indicates  possible  instrument  output  errors  which 
result  from  the  intercoupling  phenomena. 

Considerations  of  damping  constants  discussed  in  Appendix  D 
require  the  two-deqree-of-f reedom  analysis  enclosed  herein.  The  simi- 
larity between  lightly  and  highly  damped  instruments  is  emphasized 
(see  paragraph  C.7).  As  will  be  discussed  further,  the  lightly  damped 
anoroach  allows  excellent  instrument  response  about  the  three  rotational 
axes  and  enables  rapid  radial  and  axial  translation  which  permits  rapid 
warmup  while  facilitating  instrument  fill. 


C.2  SCOPE 


For  a simplified  linearized  model,  the  angular-moment  coupling  of 
the  control  loops  about  the  input  and  output  axis  of  a gyroscopic  device 
is  studied.  The  analysis  is  presented  in  a general  framework  so  that 
results  are  applicable  to  gyros,  PIGAs,  and  multisensors.  Specifically 
the  following  matrix  equation  is  studied: 


’ SI_  + C_  ! H 

r se/| 

I I j 

i 

T 

! SI0  + C0 

• 1 J 

• 

w 

CD 

o 

(C-l) . 


where: 

e - float-case  misalignment 
H ■ angular  momentum 
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I = moment  of  inertia 

C = damping  constant 

S * Laplace  transform  operator 

M = externally  applied  torque  such  as  Hu>  (gyroscopic  moment) 
or  Pa  (pendulosity) 

M_  * rebalance  torque  applied  by  suspension  (Mr  is  usually  a func- 
tion of  the  float-case  misalignment) 

Subscripts 
I - input  axis 
O = output  axis 

For  most  single-degree-of- freedom  analyses,  the  two  equations 
are  treated  separately,  i.e.,  H is  assumed  not  to  affect  the  loop 
performance.  This  report  studies  in  detail  the  effects  of  various 
control  strategies  for  M on  the  instrument  performance  described  by 
(C-l). 

C.3  JUSTIFICATION  FOR  NEGLECTING  INERTIA 

The  eigenvalues  (poles)  of  the  open-loop  system  are 


and 


For  typical  instruments 


X - 0,  0 

■ 1.6  x 10 5 dyn-cm-s 
x 10*  dyn-cm-s 

■ 200  gm-cm2 


X - -800  s-1,  -8  x 10*  s“l,  0.,  0. 


(C-2) 
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Since  gyros  and  accelerometers  are  usually  tuned  for  a closed- 
loop  bandwidth  in  the  vicinity  of  25  Hz  (157  rad/s)  as  a first-order 
approximation,  the  fast  open-loop  poles  may  be  neglected.  In  root  locus 
terms,  this  assumption  mav  be  diagrammed  as  in  Figure  C-l.  Since  the 
fast  poles  are  far  from  the  region  where  the  open-loop  poles  will  be 
moved  for  closed-loop  operation,  the  fast  poles  can  be  omitted  as  a 
first  approximation.  Omitting  the  fast  poles  (implies  neglecting  IQ 
and  Ij  in  Equation  (C-l) . 

Similarly,  consider  a typical  lightly  damped  situation. 


H ~ Cj  = 1.6  » 10 5 dyn-cm-s 

CQ  = 1.6  x lo2  dyn-cm-s 

IQ  1^  = 200  gm-cm2 

A = -400  i 693 j , 0. , 0. 

The  first  eigenvalue  has  an  undamped  natural  frequency  of  800 
rad/s;  thus,  although  the  output  axis  is  lightly  damped,  the  first 
poles  are  far  from  the  region  of  interest  (Figure  C-lb)  and  for  a first 
approximation  both  IQ  and  Ij  can  be  neglected  (although  mid-range  points 
exist  where  this  assumption  is  invalid) . Note  that  this  result  contra- 
dicts single-loop  theory  which  would  say  the  pole  for  the  output  axis 
is  Cq/Cj  - 1 rad/s. 

C.  4 REDUCED  MOTION  EQUATIONS 


Neglecting  the  inertias  in  Equation  (C-l)  and  solving  for  the 

angles 


d 

dt 


-H 


MI  ' MRI 


Mo  - mro 


C0CI  + H2 


(C-3) 


The  above  is  a fundamental  demonstration  of  the  coupling  between 
input  and  output  axes;  that  is,  torques  applied  about  both  the  input 
and  output  axes  affect  both  6j  and  0Q.  For  H small  compared  to  Cj  or 
CQ,  the  equations  for  0j  and  0Q  decouple,  a phenomenon  discussed  further. 
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X 

B.  LIGHTLY  DAMPED 


Figure  C-l.  Schematic  root  locus  justifying  omission  of  inertia 
from  motion  equation.  Poles  at  origin  will  only  be 
moved  in  shaded  area  far  from  fast  poles. 


C. 5 GENERAL  EQUATIONS  FOR  TWO-DEGREE-OF-FREEDOM  WITH  CONTROL 

Consider  the  rebalance  torques  described  by  the  control  laws 


1 

to 

M 

K 

X 

II 

'Gn  (S) 

GIO  ,S)‘ 

’»!  IS)' 

X 

» 

0 

CO 

1  

.G0I  <S> 

Goo  ,s>. 

_ e0  <s>_ 

For  generality,  each  rebalance  torque  is  a function  of  both 
angles.  Precise  expressions  for  the  G's  are  discussed  later.  Insert- 
ing (C-4)  into  (C-3),  the  transfer  matrix  relating  the  angles  to  the 
torque  is 
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Substituting  Equation  (C-5)  into  Equation  (C-4)  determines  the 
rebalance  torques  as  functions  of  the  input  torques. 


S‘C0GII*»GI0>  * GII°00  - GIOGOI 


S(-HGII«IO(V 


RI 


-s(coGoi+HGoo) 


s(ciGoo-HGoi)  + GnGoo  - GIOGOI. 


Lmroj 


DEN 


(C-6) 


where 


DEN  is  the  demoninator  of  Equation  (C-3)  . Also 
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-M 


RI 


-M 


RO 


x 


Sa(CjC0*H»)  ♦ S(ClGoo-HG01) 


-s(hgii+gioci) 


-^oSi^oo* 


S 2 (CjCq+H2 ) ♦ S(ClG00+HGI0) 


DEN 


(C-7) 


At  this  point,  several  general  observations  are  offered.  In 
general,  the  float-case-angle  misalignments  and  both  rebalance  torques 
are  influenced  by  torque  inputs  about  both  axes.  Later  discussion  will 
demonstrate  some  simplifications  which  occur  in  practical  situations 
and  will  discuss  the  implications  of  the  cross-coupling  on  instrument 
performance. 

Note,  from  Equations  (C-3)  and  (C-4)  or  (C-5) , if  GIQ  * -SH,  or 
Gqi  « SH,  the  equations  decouple.  Specifically,  the  denominator 
becomes : 


DEN  - (SCj+Gjj) (SCq+Gqq)  (C-8) 

and  the  poles  for  each  loop  can  be  determined  independently.  Also, 

Gjq  » -SH  eliminates  the  influence  of  MQ  on  0j  (and  similarly  for  Gqi, 

Mj , and  eQ) . 

Another  option  is  the  elimination  of  input  MQ  effects  on  rebalance 
Mri  by  judicious  selection  of  the  crossfeedback  compensator  GIQ,  Equa- 
tion (C-6)  or  (C-7) . in  this  situation,  the  system  polos  are  not 
decoupled  as  in  Equation  (C-8) . 

As  will  be  discussed  further,  a situation  that  appears  frequently 
in  instrument  design  is  zero  crossfeedback;  that  is  GIQ  * GQI  * 0,  for 
which 

DEN  - S*(CjC0+H*)  ♦ SICjG^Ojj)  + (C-9) 
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For  H2  small  compared  to  CjCq,  the  loops  decouple;  however,  the 
HS  terms  in  (C-l)  , (C-3) , (C-6)  and  (C-7)  indicate  that  HOj  is  acting 
like  an  input  in  the  8Q-M0  loop  (and  similarly  H§q  in  the  Oj-Mj  loop) 

C.6  MULTISENSOR  ACTIVE  CONTROLS 

A numerical  example  is  presented  to  demonstrate  the  performance 
of  an  optimized  multisensor  where  both  the  input  and  output  axes  are 
controlled  actively  with  loops  whose  bandwidth  is  approximately  22  Hz 
Both  loops  are  equipped  with  integral  compensation  so  that  under  no 
coupling  the  dynamic  response  of  the  loops  are  identical.  Parameter 
values  are  listed  in  Table  C-l. 

TABLE  C-l.  VALUES  FOR  INTEGRAL- INTEGRAL  STUDY. 


- 

1.6  x 

10 5 dyn-cm-s 

H 

BE 

n 

10s  dyn-cm-s 

j 

^ 1.6 

10 4 dyn-cm-s 

I 

i 1.6 

10 5 dyn-cm-s 

CI 

- 

1.6 

10 6 dyn-cm-s 

| 1.6 

107  dyn-cm-s  (used  in  plots) 

1.6 

10*  dyn-cm-s 

Cj (1.41 

x IOMb2  (S+22.8  s"1) 

GII 

S (S+173  s-‘) 

V1-41 

x 10  )s-2  (S+22.8  s"1) 

Goo 

S (S+173  s-1) 

Table  C-2  compares  system  eigenvalues  as  the  deunping  about  the 
input  axis  Cj  .varies.  H2/CICQ,  less  than  0.01,  is  required  for  open 
loop  design  to  give  accurate  pole  placement  in  cross-coupled  design. 
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TABLE  02.  POLES  AS  Cj  CHANGES  FOR  INTEGRAL- INTEGRAL  STUDY. 


CI 

(dyn-cm-s) 

H2 

c c 

O 

Poles 
(s- 1 ) 

oo 

0 

-34.6  -34.6 

-69.2  ±67. 3j 

-69.2  ±67. 3j 

1.6  x 10s 

0.001 

-34.4  ±0.9j 

-72.8  ±66. 7j 

-65.8  ±67. 6j 

1.6  x 107 

0.01 

1 

-80.7  ±65. 7j 

-58.4  ±68. 5j 

1.6  x 106 

KSa 

-31.0  ±6.782j 

S2K39 

-38.4  ±69. 6j 

1.6  x 105 

1.0 

-23.4  ±9 . 65 j 

-145  ±48. 6j 

-4.2  ±58. 2j 

For  H2/CjCq  = 0.01,  the  system's  response  for  steps  in  input 

torque  are  plotted  in  Figures  C-2  and  C-3.  Although  the  responses  for 
0 /M_ , 0 /Mn,  M /M  , and  M /M  are  essentially  those  of  uncoupled 
loops,  the  coupling  through  the  numerators  of  the  transfer  matrices  is 
clearly  evident.  In  particular,  the  rebalance  torque  about  the  output 
axis  caused  by  torques  about  the  input  axis  (Hu>Q  ■ M^.)  must  be  con- 
sidered in  designing  systems  which  employ  multisensors.  The  active 
suspension  has  greatly  reduced  the  permissible  angular  motion  about  the 
input  axis  and  eliminated  the  slow  time  constants  associated  with  the 
high  rotational  damping  about  the  output  axis  and  the  relatively  soft 
maqnetic  suspensions.  The  output  axis  error  is  similar  in  both 

cases. 


C.  7 LOW  DAMPING  CONFIGURATION 

With  the  inertia  terms  neglected,  the  coupled  equations  of  motion 
are  symmetric  in  that  the  roles  of  damping  and  angular  momentum  can  be 
exchanged.  Specifically,  if  H2  > 100  CjCg,  let  Gjj  - Gqq  * 0 and  use 
the  crossfeedback  compensators  GIQ  and  GQI.  The  equations  of  motion 
are  then  identical  to  those  for  the  examples  of  (C-G)  where  the  high 
damping  CjCq  > 100  H2 . 

Section  3.5  discusses  the  difficulty  of  obtaining  suitable 
dynamic  performance  for  loops  that  have  large  damping  constants.  For 
low  damping  (H2  > 100  CjCg)  the  roles  of  damping  and  angular  momentum 
have  been  exchanged  in  the  motion  equations. 
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Practically,  suitable  dynamic  control  can  now  be  designed  for  all 
three  rotational  loops.  If  a viscosity  is  selected  so  that  the  damping 
about  the  spin  axis  is  in  the  desired  range  (about  2 x 10!  dyn-cm-s) , 
the  input-output  axis  pair  will  be  in  the  lightly  damped  region  where 
the  angular  momentum, not  the  damping  constant, dominates  the  dynamic 
response.  Since  H near  2 x 10s  dyn-cm-s  is  easily  achieved  (where  size 
18  numbers  have  been  used  for  example) , desirable  dynamic  response  can 
be  achieved  about  all  three  axes.  Also,  the  low  viscosity  fluid 
(approximately  1/100  to  1/1000  of  that  commonly  used  in  single-degree 
instruments)  will  permit  large  reduction  in  the  time  associated  with 
radial  and  axial  pull-in. 
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steps  in  input  torque 


Figure  C-3.  Coupled  integral  loops.  Rebalance  torques  for  steps  in  input  torques. 


APPENDIX  D 


ACCELEROMETER  RESOLUTION  CONSIDERING  ELASTIC  RESTRAINT 
AND  SIGNAL  GENERATOR  NOISE 


D.l  SUMMARY 

The  uncertainty  in  the  measured  torques  caused  by  uncertainty  in 
the  accelerometer  signal  generator  is  investigated.  The  effects  of 
elastic  restraint  in  the  forcers  are  also  included.  The  multisensor 
feasibility  instrument  is  used  as  a baseline.  The  signal  generator  is 
modeled  as  white  noise  with  an  rms  output  of  1 milliarc  s measured  in  the 
0.1  to  50  Hz  bandwidth. 

Numerical  examples  based  on  the  feasibility  instrument  demonstrate 
that  elastic  restraint  should  not  be  a significant  problem  in  multisensor 
accelerometer  channels  provided  that  the  radial  transducer  null  does  not 
shift,  and  that  damping  about  a measurement  axis  should  be  less  than 
2 x 107  dyn-cm-s  for  an  18-size  instrument. 

For  preliminary  calculations,  an  approximate  calculation  demon- 
strates that  50  ppm  accelerometer  performance  is  equivalent  to  0.7  cm/s 
velocity  error  for  aircraft  application. 

D.2  ELASTIC  RESTRAINT 

A simple  model  of  a clapper  magnet  pair  arranged  in  push-pull 
appears  in  Figure(D-l)  Retaining  first  order  terms  (no  second  order 
terms  exist,  but  third  and  higher  occur) , the  magnet  force  in  the  up- 
ward direction  ist 


F 


)N*ApIp 


where t 

N - number  of  turns  in  each  coil 
Ip  - current  in  primary 


( D— 1 ) 


I 
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G *=  normal  air  gap 

X » displacement 

Ap  » area  of  one  pole  face 


r 

]_ 


fl  + K 


Figure  D-l.  Schematic  of  push-pull  clapper  forces. 


The  force  is  converted  to  torque  about  the  float  center  (T)  by  multi- 
plying suspension  to  center  distance.  Equation  (D-l)  may  be  rewritten 

as: 


T - kjl  - kxX  (D-2) 

For  the  original  multisensor  feasibility  instrument  with  primary 
current  of  13.1  mA  (two  gravity  accelerometer): 


kj  ■ 95.4  dyn-cm/mA 

kY  « 8.2  x 10 4 dyn-cm/cm 

By  compensating  the  accelerometer  signal  generator  primary,  the  elastic 
restraint  can  be  reduced  by  two  orders  of  magnitude;  thus  for  design 
extrapolation: 
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10s  dyn-cm/cm 


is  also  used. 


D.  3 ACCELEROMETER  EQUATIONS  OF  MOTION 

The  following  assumptions  are  made  in  considering  the  resolution: 

(1)  Translational  motion  of  the  float  is  ignored. 

(2)  Two  identical  forces  are  assumed  so  that  the  torque  exerted 
at  one  end  is  multiplied  by  two  to  obtain  the  total  torque 
exerted  by  the  magnetic  suspensions. 

(3)  Interactions  about  other  rotational  axes  are  ignored,  so 
that  only  a single  accelerometer  channel  need  be  studied. 
This  is  close  to  valid  for  torque  summations  about  the  spin 
axis. 


(4)  The  rotational  deunping,  C,  is  high  so  that  the  float's 
inertia  is  negligible. 

(5)  Uncertainty  sources  are  actual  float  motions  which  are  caused 
by  uncertain  torque  inputs  (M^)  and  electrical  noise  (xN>  in 
signal  generator. 

(6)  The  suspension  is  tuned  by  selection  of  the  feedback  compen- 
sator cH(s).  The  reasons  for  including  a factor  equal  to 
the  damping  will  be  explained  later. 

With  the  above  assumptions,  the  system  can  be  modeled  by  the  block 
diagram  of  Figure  D-2.  The  following  quantities  are  defined. 


M ■ the  applied  torque  to  be  measured 
Mjj  * noise  torques  applied  to  float 
ME  ■ torque  error  caused  by  elastic  restraint 
M * torque  applied  by  rebalance 

Mm  * measurement  torque.  Ideally  kx  = 0 so  that  measuring  the 
” current  gives  the  correct  rebalance  torque  at  steady-state. 

M_  ■ total  torque  applied  to  float 

c * rotational  damping  coefficient 

6 » float  angular  deviation  from  null 

£ - length  of  float 
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The  measured  torque  is: 


Vs) 


kjk^H (S+a) (M+Mn)  + 

2kjk 

s2  + s| 

[a  ■ 

k \ ak  l 
X \ X 

Cl  C . 

XN 

S 3 +s 2 

r k 1 1 

+ S 

• 

,kI] 

k l ‘ 
- -f  B 

+ kjk^fca 

(D-5) 


D.  4 


LOOP  COMPENSATION 


Before  proceeding  to  investigate  the  effects  of  noise  sources  and 
elastic  restraint,  the  compensator  H(S)  must  be  designed.  For  the 
multisensor  feasibility  with  kj  of  paragraph  D.2  and  l = 5.6  cm  the 
following  terms  complete  the  torque  loop: 

k?  = 26.4  mA/cm-s- (dyn-cm-s) 

a - 22.8  s'1 

B « 173  s“‘ 

With  zero  elastic  restraint,  the  Bode  plot  for  the  closed  loop 
accelerometer  appears  in  Fig.  D-3. 

NOTE:  The  damping  constant  C has  cancelled.  (The  physical 

C in  the  forward  path  is  cancelled  by  the  compensator's 
mathematical  C.) 

The  system  has  a bandwidth  of  approximately  22  Hz  (as  described  in 
Cnapter  2,  accelerometer  loop  bandwidth  was  reduced  to  0.8  Hz). 

D.5  DISCUSSION  OF  TRANSFER  FUNCTIONS 

For  systems  with  typical  bandwidth  and  damping  constant,  the 
elastic  restraint,  kx»  has  only  a small  effect  on  the  poles.  Assume 
C * 2.5  x io 5 dyn-cm-s  and  parameters  used  in  feasibility  instrument. 
Referring  to  (D-5)  and  performing  the  calculations: 


k l 

B - -S- 


173  - 2 ~ 173  ~ B 


kiV 


k *B 

x 


1.41  x io*1  - 320  ~ 1.41  x io1*  ~ kjk£f 
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Thus,  for  systems  with  higher  damping,  elastic  restraint  can  be  ignored 
in  the  denominator.  (Accelerometer  loops  will  have  higher  damping  con- 
stants.) If  kx  affects  poles,  a,  8,  and  ki  can  be  adjusted  to  maintain 
pole  position.  Thus,  for  systems  where  inertia  can  be  neglected,  if  a 
desired  response  is  given  by  H(S)  and  the  physical  C is  changed,  one 
only  need  change  the  constant  in  the  feedback  loop  proportionately. 

For  a given  pole  placement  (i.e.,  given  HC(S)),  the  float  position 
is  proportional  to  the  applied  torques  but  inversely  proportional  to 
the  damping  constant.  High  damping  results  in  small  motion  so  that  the 
motions  to  be  measured  become  comparable  to  the  SG  noise  and  resolution 
problems  must  be  considered.  The  float's  actual  position  is  propor- 
tional to  the  SG  noise  input  with  no  dependence  on  damping.  The  float's 
position  does  not  depend  on  elastic  restraint. 

Referring  to  the  measured  torque  (D-5) , the  measured  torque  is 
proportional  to  the  applied  torques  with  no  dependence  on  damping  con- 
stant. The  interesting  term  is  the  dependence  of  measured  torque  on 
the  signal  generator  noise  input.  With  typical  elastic  restraint,  the 
measured  torque  caused  by  SG  uncertainty  is  proportional  to  the  damping. 
For  typical  systems 


and  higher  derivatives  have  a larger  effect  on  rebalance  noise  (see 
paragraph  D.5). 

Since  proposed  accelerometer  loops  could  have  rotational  damping 
three  to  four  orders  of  magnitude  higher  than  conventional  measurement 
loops  unless  consideration  is  given  to  float  shaping,  the  influence  of 
the  signal  generator  noise  and  the  elastic  restraint  will  be  studied 
further. 

0. 6 MEASURED  TORQUE  UNCERTAINTY  CAUSED  BY  SG  NOISE  IN  THE  FEASIBILITY 

INSTRUMENT 

Consider  the  torque  summation  about  the  spin  axis  which  is 
essentially  decoupled  from  motions  about  the  other  axes.  The  following 
values  are  summarized  for  the  multisensor  feasibility  instruments 
originally  designed. 
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a * 22  s“l 

6 * 173  s"1 

kx  = 8.2  x 10“  dyn-cm/cm 

kj  * 95.4  dyn-cm/mA 

j.  — r a n\A  

Kl  * cm-s  (dyn-cm-s) 

C ■ 2.5  * 10  s dyn-cm-s 

SG  noise  with  no  float  motion  was  experimentally  described  in 
measurements  from  0.1  Hz  to  50  Hz  as  white  noise  with  spectrum  of  the 
float  angle 


♦0(S)  * 3.3  x 10-“  rad2/rad/s  (D-6a) 

which  corresponds  to  1 milliarc  s measured  from  0.1  to  50  Hz.  convert  to 
position  from  angle  with  f - 5.6  cm  gives  the  spectrum  of  xN 

*x<S)-  2.59  10-“  jgjy  (D-6b, 

With  known  input  spectrum  and  transfer  function,  the  output  spec- 
trum can  be  determined.  From  the  output  spectrum,  rms  values  of  the 
expected  values  are  determined.  Table  D-l  tabulates  the  rms  measurement 
torque  uncertainties  for  several  conditions . Deunping  ratios  of 
C * 2.5  x 10s  and  2.5  * 10*  dyn-cm-s  are  compared  for  models  with  the 
estimated  elastic  restraint  emd  zero  elastic  restraint. 

The  numerical  examples  which  were  tabulated  for  a loop  whose  band- 
width is  22  Hz  demonstrates  that  the  rms  rebalance  noise  caused  by  the 
signal  generator  uncertainty  is  proportional  to  the  deunping  constant 
for  a given  bandwidth.  The  presented  data  can  be  transformed  to  other 
bandwidths  since  the  rms  torque  noise  is  proportional  to  the  bandwidth 
to  the  three  halves  power.  For  deunping  constants  of  2.5  * 10*  dyn-cm-s 
which  is  representative  of  the  feasibility  instrument,  the  elastic 
restraint  has  little  effect  on  the  rebalance  torque  noise.  For  the 


presently  configured  feasibility  instrument  with  0.8  gm-cm  pendulosity, 
the  rms  noise  is  0.12  gravity,  a relatively  high  number  which  leads  to 
the  following  questions: 

(1)  What  are  the  effects  of  the  torque  noise  on  stability  testing 
of  the  feasibility  instrument? 

(2)  What  are  the  implications  for  actual  instrument  design? 


TABLE  D-l.  RMS  TORQUE  UNCERTAINTIES  CAUSED  BY  SG  NOISE. 
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Because  of  the  coupling  between  the  gyro  torquer  and  the  radial 
transducer,  the  effective  noise  in  the  signal  was  larger  than  that  used 
in  the  above  analyses.  The  loop  bandwidth  was  reduced  to  0.7  Hz  and 
the  loop  output  was  filtered  externally  by  an  analog  filter  with  0.02  Hz 
bandwidth. 

Very  low  frequency  or  steady-state  null  shifts  do  not  appear  in 
the  spectra  (D-6a) . Using  (D-5)  to  relate  the  measured  torque  to  the 
signal  generator  null  output  at  steady-state: 

MM  " 2kXXN  (D"7) 

With  a pendulosity  of  0.8  gm-cm  and  a design  goal  of  100  pg  bias 
stability,  the  drift  in  the  measured  torque  must  be  less  than  0.08 
dyn-cm.  To  allow  for  other  bias  sources,  the  torques  calculated  from 
(D-7)  should  be  less  than  0.08.  For  this  discussion  0.05  dyn-cm  is 
selected.  With  kx  * 8.2  x 104  dyn-cm/cm  (the  feasibility  instrument 
tuned  for  2g  with  no  signal  generator  tuning),  signal  generator  drift 
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of  less  than  3 x 10" 7 cm  is  required.  With  k - 10 5 dyn-cm/cm  (the 
feasibility  instrument  with  signal  generator  tuning),  3 x io_s  cm  is 
required.  Note  that  these  results  have  no  dependence  on  the  damping 
constants.  With  series  parallel  tuning  (SECTION  2),  null  shift  (AX„) 
can  be  related  to  capacitor  changes  (AC)  by  the  approximate  relation 


where 

CQ  is  the  nominal  capacitance 

gQ  is  the  nominal  magnet  pole  face  clearance 

For  the  feasibility  instrument  g,  is  0.015  cm  and  C is  approximately 
0.038  pF.  Thus  capacitance  changes  must  be  less  than  0.76  pF  and  76  pF 
for  the  untuned  « id  tuned  signal  generator  respectively.  These  capaci- 
tor changes  are  relative  changes.  If  all  capacitors  change  identically, 
no  shift  occurs.  The  implications  of  capacitor  changes  are  discussed 
in  paragraph  3.6. 

The  following  considerations  of  the  signal  generator  noise  on 
multisensors  for  flight  use  are  discussed. 

(1)  Elimination  of  capacitance  sensitivity  (this  has  been 
covered  in  the  previous  discussion  of  feasibility  tests) . 

(2)  Effects  of  elastic  restraint  on  velocity  error;  and 

(3)  The  trade  between  damping  constant  and  loop  bandwidth. 

Since  accelerometers  are  used  to  determine  velocity,  accelerometer 
specifications  should  be  converted  to  velocity  so  that  meaningful 
evaluation  may  be  accomplished.  If  one  desires  50  ppm  accelerometer 
response  and  if  a nominal  operating  velocity  of  300  mile  per  hour  is 
assumed,  the  velocity  error  must  be  less  than  0.7  cm/s  (50  ppm  multi- 
plied by  300  mile  per  hour) . 

Velocity  errors  caused  by  elastic  restraint  are  discussed. 
Referring  to  Figure  D-2  and  Equation  (D-4)  and  (D-5) , and  omitting 
signal  generator  and  input  noise;  the  error  in  the  measurement  torque 
is  derived: 
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Me  = M - Mm 


(D-8) 


kvae 


ky*e  i 

- -g—  J S + kjk^a 


/ 

For  steady-state  torque  input  (acceleration  input)  , the  measurement 
error  is  zero  at  steady-state;  however,  the  integral  of  torque  is  pro- 
portional to  velocity.  For  a step  input  in  M,  (corresponding  to  one 
gravity)  the  steady-state  velocity  error  is: 


B k (P.) 

Steady-state  velocity  error  = A— p- — 

01  K jK 

where  the  torque  M has  been  written  as  pendulosity  times  acceleration. 
With  the  tuning  for  22  Hz  bandwidth  and  compensator  elastic  restraint 
of  1000  dyn-cm,  the  steady-state  velocity  error  for  a one-gravity  input 
in  acceleration  is 


Velocity  steady-state  - 10 1 (dyn-cm-.) 

with  damping  constant  of  5 * 10 5 dyn-cm-s  or  greater,  the  steady-state 
velocity  error  is  less  than  0.005  cm/s  which  is  much  less  than  the 
0.7  cm/s  equivalent  to  50  ppm  navigation. 

The  rebalance  torque  noise  resulting  from  the  signal  generator 
noise  is  estimated  at  94  dyn-cm  (Table  D-l)  or  0.12  g with  0.8  gm-cm 
pendulosity.  This  calculation  is  based  on  the  22  Hz  measurement  band- 
width and  the  signal  generator  spectrum  of  (D-6) . For  the  given  loop 
dynamics,  the  integral  of  the  rebalance  torque  is  proportional  to  the 
flight  velocity.  With  0.8  gm-cm,  the  rms  velocity  is  1.1  cm/s  which 
is  greater  than  the  0.7  cm/s  comparable  to  50  ppm  performance.  In 
addition,  calibration  is  difficult  because  the  rms  noise  is  three 
orders  of  magnitude  greater  than  the  calibration  accuracy  (50  ppm). 

Since  bandwidth  reduction  is  undesirable  for  operational  systems,  the 
multisensor  transverse  axes  should  be  designed  with  damping  an  order 
of  magnitude  smaller,  a suggestion  which  is  explored  in  paragraph  3.8  and 
Section  4. 
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APPENDIX  E 


CONVECTION  TORQUES  ABOUT  TRANSVERSE  AXIS 


E.l  SUMMARY 

This  Appendix  demonstrates  the  convection  pattern  resulting  from 
an  end-to-end  thermal  gradient  and  the  effects  on  accelerometer  stability 
in  the  multisensor.  For  the  multisensor  feasibility  instruments  ther- 
mally induced  torques  of  the  order  of  5 dyn-cm/g-°F  (°F  of  end-to-  o- 
end  thermal  difference)  can  be  expected. 

E.2  DISCUSSION 

An  approximation  of  the  effects  of  an  end-to-end  thermal  gradient 
in  the  multisensor  is  obtained  for  the  model  of  a cylindrical  float 
depicted  in  Figure  E-l.  Similar  analyses  have  been  conducted  about  the 
cylinder  axis.  The  following  assumptions  are  made: 

(1)  Fluid  inertia  may  be  neglected; 

(2)  Using  the  coordinates  of  Figure  E-l,  a linear  thermal  gradient 
exists  along  the  float  so  that  the  temperature  T(Z)  at  any 
point  in  the  fluid  is  given  by: 

T = T(Z)  = TQ  (1  - 0|)  (E-l) 

where : 

Tq  « mean  temperature; 

AT  * temperature  difference  between  ends 

(3)  The  cylindrical  float  is  coaxial  with  the  cylindrical  case; 

(4)  Gravity  is  perpendicular  to  the  float's  axis  as  shown  in 
Figure  E-l; 
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Figure  E-l.  Coordinates  for  convection  analysis. 

(5)  The  fluid  flow  can  be  divided  into  three  regions  of  interest: 
region  1 is  an  annulus  of  radius  R,  length  l and  thickness 
h.  Flow  in  the  annulus  is  in  the  Z direction  only,  an  assump- 
tion whose  validity  will  be  demonstrated  later.  Assuming  h 
much  smaller  than  l and  R,  the  pressure  gradient  in  the 
annulus  is: 


where: 

P 

Q 

u 


(E-2) 

j 


1 

j 

pressure  is  a function  of  axial  position  Z and 
circumferential  angle  9 


volume  flow/unit  length  in  Z direction  and  a 
function  of  6 


absolute  viscosity  which  is  assumed  uniform 
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(6)  The  end  volumes  are  large  so  that  viscous  effects  are  small. 
The  pressure  gradient  then  depends  upon  gravity.  With  the 
temperature  gradient  defined  by  (E-l),  the  pressure  gradient 
in  region  2 is: 

H ■ p09  (3  * "T”) 

where: 

g = specific  force 
k = fluid  volume  gradient 
C0  = fluid  density  at  reference  condition 
In  region  3,  the  pressure  gradient  is  given  by: 

1st  • V i1  - ¥-)  (E-4> 

Proceed  to  solve  for  the  fluid  pressure.  Circumnavigating 
the  float  by  the  path  depicted  in  Figure  E-l,  the  path 
integral  must  be  zero: 


P div  P • dx  = 0 (E-5) 

Substituting  (E-2) , (E-3)  and  (E-4)  into  (E-5) , allows  solution  for 
Q<e> 


h3  k AT  Png  R cos  0 

0 <e>  (E'6) 

Integration  of  (E-4)  gives  the  pressure  gradient,  at  the  left  end: 

l ( VT  \ 

p (-  f,  0»  R)  = pq  + pQg  R ^l  + — J (l  - cos  e)  (e-7) 


where 


X = R cos  0 

and  where  PQ  is  a reference  pressure  which  need  not  be  determined  for 
this  work. 
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Integration  of  (E-2)  with  the  boundary  condition  (E-7)  gives  the 
pressure  in  the  annulus: 


P(0,z) 


p0g 


1 + 


kAT 

2 


- pog 


R COB 


r VT  1 

+ L"r  pogR  COB  eJ! 


(E-3) 


Note  that  (E-8)  is  a steady-state  solution  for  the  pressure 
distribution  and  that  (E-6)  is  convection  path.  Return  and  verify 
assumption  5 where  circumferential  flow  in  the  annulus  was  assumed 
zero.  For  flow  along  the  circumference: 


12  pQr 
h1 


1 3P 
R 96 


k AT  \ 

— z)  g sin  0 


(E-9) 


Substitution  of  (E-9)  into  (E-8)  gives  0R  equal  zero  which  verifies  the 
assumption. 

Since  the  pressure  distribution  and  the  flow  are  known,  the 
forces  and  moments  acting  upon  the  float  can  be  determined. 

The  force  in  the  vertical  (X)  direction: 


F * 

rx 

pQg  it  R2  i 

(8.10) 

In  the  (Y)  direction 

F m 

X 

0 

(E.ll) 

and  in  the  axial  direction 

I 

1 

1 

M 

1 

0 

(E . 12) 

Since  the  flow  for  positive  X is  the  negative  of  the  flow  for  negative 
X (the  cos  8 in  (E-6) , the  solution  for  the  bottom  half  is  antisymmetric 
with  respect  to  the  solution  for  the  upper  portion  so  that  (E-12)  results. 

From  the  pressure,  the  torque  induced  by  normal  pressure  forces  is 
derived.  For  this  symmetric  example,  all  torques  are  about  the  Y axis. 

The  contribution  from  the  ends  is: 
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( E— 1 3 ) 


Mend  _ * kT  pOg 

AT  4 


The  torques  generated  by  the  normal  pressures  in  the  annulus  is : 
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ann 

AT 


12  POg  R2 


l2  K 


(E-14) 


The  torques  caused  by  the  pressures  in  the  annulus  and  in  the  ends  are 
in  opposite  directions,  thereby,  tending  to  cancel  one  another.  Since 
the  ends  were  assumed  large  so  that  flow  rates  were  small,  no  viscous 
shear  force  is  contributed  in  the  end  regions.  In  the  annulus,  the 
shear  stress  exerted  by  the  fluid  on  the  float  is: 


T 


z 


6u  Q (6) 
h2 


and  the  torque  is: 
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(E-15) 


(E-16) 


Estimate  the  effect  of  thermal  gradient  on  a multisensor  accelerometer 
where  the  following  parameters  are  from  the  feasibility  instrument: 


R - 2 cm 

t ■ 4 cm 

kT  - 4 k 10- 

g • 980  cm/s2 

pQ  ■ 2.655  gm/cm* 

The  estimated  torques  are: 

^ - 13  dyn-cm/°F 


Mann 

At 


17  dyn-cm/°F 


147 


The  viscous  torques  are  negligible  since  h<*£  or  R.  With  h = 0.0254  cm 
the  viscous  torque  is: 


M 


vis 

AT 


0 . 3 dyn-cm/°F 


With  pendulosity  of  800  dyn-cm  at  one  gravity,  the  scale  factor  sensi- 
tivity is: 


SF 


5000  EEH 

°F 


APPENDIX  F 


STATIC  CROSS-COUPLING  IN  MULTISENSORS 


F.l  SUMMARY 

Based  on  physical  phenomena  and  misalignment  of  forces  and 
torques,  an  error  model  for  multisensor  1 g and  earth-rate  tests  is 
formulated.  The  model  will  be  used  for  locating  instrument  axes  and 
determining  coupling  between  accelerometer  and  gyro  functions. 

The  model  indicates  that  large  pendulous  unbalance  along  OA 
appears  as  acceleration-sensitive  drift  which  previous  multisensor  tests 
demonstrated  to  be  in  the  range  of  typical  gyros. 

For  calibrating  the  accelerometer,  precise  knowledge  of  the 
current-torque  relation  is  not  required;  that  is,  signal-acceleration 
calibration  can  be  performed  directly.  The  accelerometer  output  con- 
tains a scale  factor  which  should  be  set  to  tune  out  translational 
motion  effects.  Cross-coupling  effects  will  determine  the  position  of 
the  accelerometer  IA  so  that  the  axis  will  not  coincide  with  IRA  or  the 
uyro  IRA. 

F.2  MULTISENSOR  MODEL  AT  STEADY- STATE 

As  a first  approximation,  only  steady-state  will  be  considered. 
Second-order  terms  will  be  neglected.  At  steady-state,  Newton's  laws 
must  be  satisfied  for  the  rotational  motion 

EM  * M..  + M + M + M_  + M_  = u>  * (F-l) 

— M — g — g 2 — R — B -S 

and  translational  motion 

EF  » 5^  + Tb  + £r  “ 0 (F-2) 
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where: 

F_  * buoyancy 

— D 

Fw  = magnetic  suspension  force 
F ■ reaction  force 

R 

M_  ■ buoyant  torque 

Mg  * acceleration-proportional  torques 

M = acceleration-squared  torques 

g2 

“m  * magnetic-suspension  torques 
Md  * reaction  torque 

J\ 

a:  * angular  velocity  of  float  qj/IS 

Hg  ■ angular  momentum  of  wheel 

Note,  since  u>  * u>E  * 7.27  * 10" ! rad/s,  many  terms  have  been 

neglected. 

Determine  the  vectors  in  Equation  (F-l)  and  (F-2)  for  the  spin, 
input  and  output  axes.  Define  the  coordinate  triad  (S,  I,  O)  . The 
rate  of  change  of  angular  motion  is  determined: 

0 J -a) 

h — 

“o|  ° 

r — 

— (*)_  | u 

1 I 

where  H is  written  with  3 components  since  float  and  reference  axes  do 
not  coincide.  HI  and  HQ  are  assumed  small. 

In  float  coordinates,  torques  proportional  to  acceleration  are 
determined  from 
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where 


a * acceleration 

Ug  ■ mass  unbalance  (small)  along  input  and  spin  axes 
P - pendulosity  (large  mass  unbalance)  along  output  axis 


The  g2  terms  are 
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Compliance  terms  are  usually  small  (15  meru/g2  for  example);  thus, 
compliance  multiplied  by  another  small  quantity  is  negligible. 

Define  the  vector  of  reaction  torques 
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and  a vector  of  buoyancy  torques 
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( F—  7 ) 


The  reaction  forces  which  are  caused  by  magnetic  misalignments  and  flex 
leads  are: 


RS 


F 


RI 


(F-3) 


RO 
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and  the  buoyancy  forces; 


-B 


BSaS 

BI3I 

Boao 


(F-9) 


where  B is  the  buoyancy  along  the  indicated  axis.  As  discussed  in 
Appendix  A,  this  representation  of  buoyancy  is  probably  an  oversimplifi- 
cation. 

Neglecting  elastic  restraints,  a valid  assumption  with  null  seek- 
ing suspensions  and  nonlinearities,  the  magnitude  of  the  forces  and 
torques  is  proportional  to  the  secondary  current  in  the  device.  Because 
of  forcer  misalignments  magnetic  anomalies,  and  float-case  misalignments, 
the  magnetic  forces  do  not  align  themselves  perfectly  with  the  float's 
S-I-O- triad.  Thus,  the  forces  and  torques  may  be  expressed  by  the 
matrix  relation: 
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(F-10) 


where  the  overline  indicates  a small  quanity  (if  forcers  were  ideal 
these  quantities  would  be  zero)  and  i - current  in  forcer  axis  indica- 
ted in  Figure  F-l. 


Equation  (F.10)  is  very  general  and  includes  all  possible  linear  cross- 
couplings. 

Incorporate  Equation  (F.l)  through  (F.10)  to  determine  the 
currents  in  the  radial  forcers  and  gyro  torquer. 
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where  the  overline  indicates  small  quantities  and 


a 


For  an  ideal  instrument,  the  overlined  terms  are  small  and  the  following 
relations  are  based  on  geometric  considerations. 


q12  - 

-q22 

(F-12a) 

q13  * 

q23 

(F-12b) 

q14  - 

"q41 

(F-12c) 

q34  " 

q44 

(F-12d) 

The  approximate  equal  signs  of  Equations  (F-12a)  and  (F-12c)  are 
caused  by  the  instrument's  pendulosity  which  causes  the  geometric  float 
center  to  differ  from  the  center  of  gravity  so  that  slightly  different 
moment  arms  exist  for  the  two  forcers. 
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Figure  F-l.  Schematic  of  forces  and  torques  in  multisensor. 


Previous  multisensor  analysis  has  assumed  that  the  accelerometer 
readout  for  acceleration  along  the  input  axis  was  the  rebalance  torque 
about  the  spin  axis  which  is  roughly  determined  by  the  radial  suspension 
forces  exerted  along  the  input  axis.  Therefore,  define  the  accelerometer 
output  signal  Sj.  as  follows 


i - k i 
1I1  *1  XI2 


(F-l 3 ) 


where  kj  is  a constant  whose  selection  will  be  discussed  further. 
Substitution  of  Equation  (F-ll)  into  (F-13)  results  in 
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where  second-order  terms  have  been  neglected. 

Expressions  similar  to  Equation  (F-14)  can  be  derived  for  the 
currents  igl  and  ig2  and  for  the  gyro  torquer  current 
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F.  3 COMMENTS 

The  position  of  the  axes  and  the  cross-couplings  at  steady-state 
do  not  depend  on  the  signal  generators  (neglecting  SG  null  shifts,  see 
Appendix  C) . At  steady-state  Newton's  laws,  Equations  (F-l)  and  (F-2) 
must  be  satisfied.  Thus,  in  Equations  (F-14)  and  (F-13),  one  is  con- 
cerned only  with  the  steady-state  current  and  not  with  the  path  the 
current  followed  to  reach  said  level. 

In  the  gyro  torquer  current  Equation  (F-15) , the  cross-coupling 
q^2  will  cause  the  instrument's  measured  input  axes  to  be  located  at  a 
position  other  than  perpendicular  to  the  spin  axis  and  to  the  axis  of 
symmetry  of  the  float.  (This  observation  pertains  to  all  gyros.) 

Because  of  the  cross-coupling  qgl  and  q52,  the  pendulosity  along  the 
output  axis  appears  as  a gyro  acceleration  drift  coefficient.  Equation 
(E-15) . Measurements  of  the  multisensor  feasibility  instrument's 
acceleration  drifts  were  in  the  same  range  as  those  for  similar  instru- 
ments so  that  one  may  conclude  these  couplings  are  small. 

Consider  the  accelerometer  signal  Sj,  Equations  (F-13)  and  (F-14). 
Any  value  of  kj  (except  that  which  zeroes  or  makes  small  the  first 
expression  in  Equation  (F-14)  will  give  acceleration  information  at 
steady-state.  Thus,  the  precise  relation  between  current  and  torque 
need  not  be  determined.  The  relation  between  current  and  acceleration 
is  sufficient.  The  current  weighting  factor  should  be  selected  so  that 
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( F— 16) 


(£*34  * kI  *44>  = 0 

This  minimizes  the  effects  of  translational  motion  on  the  accelerometer 
readout  (F-14) , while  Equation  <F-12c)  indicates  acceptable  gain  for 
Paj.  Thus,  error  torques  associated  with  a change  in  flotation  tempera- 
ture are  greatly  reduced  since  B^a^  is  usually  large  compared  to  other 
thermal  error  sources.  The  modes  associated  with  the  transient  motion 
might  be  slower  than  the  rotational  modes  (see  Appendix  B)  . Selecting 
kj.  according  to  Equation  (F.16)  eliminates  the  slow  modes  from  the 
accelerometer  signal  Sj  although  the  slow  mode  remains  in  the  float's 
actual  behavior. 

Proper  kI  could  be  chosen  byi 

(1)  Raising  the  instrument  temperature  uniformly,  and/or 

(2)  Moving  the  float  translationally  and  allowing  return  to 
equilibrium.  kj.  is  adjusted  until  the  translational  effects 
are  not  seen  in  Sj. 

The  aQ  and  ag  terms  reflect  that  the  multisensor  input  axis  will 
not  be  on  IRA  and  probably  will  not  coincide  with  the  gyro  input  axis. 
Additionally,  the  accelerometer  signals  will  contain  rate  information 
because  of  the  coupling  phenomena.  Multisensor  cross-coupling  tests 
will  attempt  to  evaluate  these  effects. 
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